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FOREWORD 


This  report  was  prepared  by  Electro-Optical  Systems,  Inc,  Pasadena, 
California  on  Air  Force  Contract  AF  33(616) -5977  under  Task  No.  7116-02, 
,lRe«.tjarch  in  Energy  Conversion  Techniques",  of  Project  No.  7116,  ’’Energy 
Con^^ersion  Research”.  The  contract  was  monitored  by  Messrs.  Ken  Vickers 
and  TE.  D.  Stephens,  of  the  Aeronautical  Research  Laboratories. 

The  investigation  began  6  June  1958  and  was  accomplished  by  the 
Adv££  xiced  Power  Systems  Division  of  EOS. 

In  its  first  phase,  as  reported  in  WADC  TN59-184,  the  work  was  done 
by  MEx,  Donald  H.  McClelland,  Project  Supervisor  under  the  direction  of 
Dr.  H.  Fisher,  Division  Manager.  In  its  second  phase,  as  reported 

herfe  ,  the  work  was  done  by  Dr.  Frederic  E.  Fuller,  Project  Supervisor 
undfcx:  the  direction  of  Mr.  Donald  H,  McClelland,  Manager,  Energy  Con¬ 
vert  ion  and  Regulation  Department.  This  is  the  final  report  of  the 
secoxid  phase  of  this  investigation. 

The  research  reported  herein  is  exploratory.  The  concept  treated 
was  drirst  proposed  by  Mr.  Elmer  Johnson  of  ARL  and  was  reported  in 
WADC  TR  56-182  by  .  Robert  Hunter. 
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ABSTRACT 


The  feasibility  of  a  solar  powered  aircraft  is  investigated. 
The  configuration  reached  is  a  high  aspect  ratio  flying  wing  with 
ramjet  thrust.  A  statement  is  made  with  respect  to  feasibility 
with  a  recommendation  of  a  program  for  further  werk.  The  statement 
of  feasibility  is  supported  by  experimental  investigation  of  solar 
radiation  collection  and  inlet  air  diffusion  and  finally  by  design 
and  performance  analyses. 
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1 .  INTRODUCTION 

The  use  of  solar  energy  as  a  power  source  has  always  been 
intriguing  because  of  its  inexhaustibility.  Since  the  advent  of 
the  space  flight  concept,  scientists  and  engineers  have  sought  methods 
for  adapting  the  solar  source  of  radiant  energy  to  the  propulsion  of 
spacecraft.  In  recent  years,  another  concept  has  gained  consideration  - 
the  solar  powered  aircraft. 

The  application  of  solar  energy  to  aircraft  propulsion  is  particularly 
interesting,  because  it  would  make  possible  a  device  that  could  circle 
the  earth  continuously,  in  the  same  manner  as  a  satellite,  only  at  a 
much  lower  altitude  and  velocity.  The  fundamental  object  of  research 
and  development  conducted  under  this  contract  was  to  investigate  problems 
concerning  the  feasibility  of  flight,  employing  a  solar  powered  ramjet. 

This  report  is  in  three  major  parts;  the  experimental  investi¬ 
gation  of  radiation  trapping,  the  experimental  investigation  of  the 
subsonic  part  of  the  inlet  air  diffusion,  and  finally,  the  design 
and  performance  analysis  development.  These  investigations  are  brought 
to  culmination  in  a  statement  on  the  technical  feasibility  of  the  solar 
powered  aircraft.  A  recommendation  is  offered  for  the  direction  of 
further  investigation. 

The  experimental  investigation  of  radiation  trapping  is  reported 
in  Section  13  of  Part  I.  The  material  of  this  section  has  been  directly 
transferred  from  WADC  TN59-184  for  ready  reference.  The  experimental 
investigation  of  the  inlet  diffuser  is  reported  in  Part  II.  In  Sec¬ 
tions  2  through  13  of  Part  I,  is  presented  the  solar  powered  aircraft 
design  development  and  performance  analysis  with  configurations  and 
estimated  performance  all  in  support  of  the  feasibility  statement  and 
the  recommendation  for  further  investigation. 

Manuscript  released  by  author  14  April  1961  for  publication  as  an 
ARL  Technical  Report 
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2.  THE  AIRCRAFT  CONFIGURATION 

The  configuration  of  the  solar  powered  aircraft  as  developed  in  this 
study  is  shown  on  Figure  1.  This  particular  configuration  emerged  in  the 
course  of  the  work  on  the  basis  of  the  following  cone idarat ion* : 

1.  High  wing  span  to  chord  length  aspect  ratio. 

2.  High  Reynolds  number  based  on  flight  speed,  the  wing  chord 
length,  and  the  free  stream  air  viscosity. 

3.  High  ratio  of  area  of  solar  beam  cross-section  to  wing  area. 

4.  Optimum  areas  of  the  internal  stream  and  of  the  free  stream 
entering  the  internal  stream. 

5.  Optimum  compromise  of  pressure  loss  in  the  internal  stream 
with  lift  and  pressure  drag  in  the  external  stream. 

6.  The  two-shock,  unsymmetrical  supersonic  diffuser. 

7.  The  diagonal  array  of  heat  transfer  elements  in  the  internal 
stream. 

8.  The  spectrally  selective  surfaces:  upper  wing  surfaces,  both 

inner  and  outer;  the  lower  wing  surface,  both  inner  and  outer  and 
the  heat  transfer  foils. 

The  high  aspect  ratio  and  Reynolds  number  tend  toward  low  induced 
and  friction  drags,  respectively.  The  high  ratio  of  solar  radiation 
beam  cross-section  to  wing  area  tends  toward  maximum  power  loading. 

The  use  of  the  entire  span  of  the  wing  leading  edge  for  inlet  to  the 
internal  stream  yields  maximum  internal  airflow  consistent  with  the  drag 
of  the  external  stream.  As  will  be  discussed  later,  this  choice  may  not 
be  optimum  with  respect  to  pressure  loss  in  the  internal  stream.  By  use 
of  the  maximum  cross-section  of  the  internal  stream  otherwise  consistent 
with  the  external  drag  and  the  internal  pressure  loss,  the  pressure  loss 
due  to  heat  transfer  is  minimized.  In  Sections  10.1  and  10.2  analyses  are 
presented  in  support  of  the  last  two  points:  the  first  of  these,  based 
on  a  consideration  of  propulsive  efficiency,  suggests  including  the  percent 
of  span  to  be  used  as  diffuser  inlet  width  as  a  design  parameter;  the  second 
of  these,  based  on  a  consideration  of  the  Reynolds  analogy,  confirms 
the  choice  of  maximum  internal  stream  cross-section. 

With  respect  to  the  compromise  of  pressure  losses  in  the  internal 
stream  with  lift  and  drag  in  the  external  stream:  these  effects  are 


o 


critically  dependent  on  the  length  and  slopes  of  the  upper  and 
lower  wing  surfaces.  In  the  design  analysis,  the  lengths  and  slopes 
of  the  inlet  diffuser,  the  heat  transfer  and  jet  nozzle  section  are 
determined  for  maximum  wing  loading  of  the  aircraft. 

The  advantages  of  the  two-shock  unsymmetrical  supersonic  diffuser 
are  discussed  in  Section  9.  These  advantages  are:  a  high  pressure 
recovery  in  the  supersonic  diffuser  and  the  delivery  of  the  inlet 
airstream  to  the  subsonic  diffuser  entrance  with  minimum  disturbance 
of  the  internal  and  external  streams. 

The  diagonal  array  of  heat  transfer  elements  in  the  internal  stream 
permits  intercepting  the  solar  beam  by  equal  or  greater  area  of  absorbing 
surface  but  with  both  foil  lengths  and  diagonal  spacing  available  to  be 
determined  for  best  design.  Of  course,  as  the  diagonal  spacing  varies 
the  slope  of  the  diagonals  also  varies  accordingly. 

The  spectrally  selective  surfaces  serve  to  trap  solar  radiation 
energy  within  the  aircraft  wing  where  it  is  delivered  to  the  internal 
stream.  The  spectrally  selective  surfaces  are  used  as  follows: 

1.  A  "green  house  effect"  on  the  wing  upper  surface  where  the 
wing  upper  surface  has  high  transmittance  for  solar  radiation 
and  high  reflectance  for  the  characteristic  thermal  radiation 
of  the  heat  transfer  foils. 

2.  Heat  transfer  foils  with  high  absorptivity  for  the  solar 
radiation  and  low  emissivity  for  their  characteristic  thermal 
radiation . 

3.  Wing  inner  surfaces  highly  reflective  to  the  heat  transfer  foil 
characteristic  thermal  radiation. 

4.  Wing  outer  surfaces  of  low  emissivity  for  their  own  characteristic 
thermal  radiation. 
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3.  THE  AIRCRAFT  FLIGHT  CONDITIONS 

The  flight  conditions  most  favorable  for  feasibility  were 
chosen.  These  conditions  are  determined  by  the  requirements: 
first,  that  the  sun  be  directly  overhead  the  aircraft,  and 
second,  that  the  altitude  be  optimum. 

Conditions  of  a  24  hour  period  earth  satellite  at  constant 
altitude  are  given  in  Table  I,  The  atmospheric  properties  used 
are  given  in  Ref.  1,  the  ARDC  Model  Atmosphere.  In  Table  I,  the 
quantities  of  special  interest  are: 

Z  -  altitude 

M  *  aircraft  flight  Mach  number 

v/V^  -  aircraft  flight  Reynolds  number  per  foot. 

/ v  \  ■  Reynolds  number  per  foot  of  the  internal  flow 

^  ^  s  2  at  the  diffuser  outlet. 

The  aircraft  optimum  performance  and  configuration  depend 
critically  on  the  flight  altitude. 


4.  INTRODUCTORY  DISCUSSION  OF  THE  AIRCRAFT  PERFORMANCE 

The  aircraft  configuration  is  shown  in  Fig.  1  and  a  summary  discussion 
of  the  basis  for  the  choice  of  this  configuration  is.  given  in  Section  2. 

In  the  following,  an  elementary  discussion  of  the  aircraft  flight  processes 
will  be  given. 

The  vector  sum  of  all  the  surface  normal  and  tangential  forces  over 
the  exterior  and  interior  surface  of  the  aircraft  yield  the  usual  equations 
for  the  equilibrium  of  lift,  total  weight,  drag  and  net  thrust.  The  net 
thrust  is,  by  definition,  the  vector  sum  of  all  aerodynamic  forces  on  the 
interior  surfaces  of  the  aircraft;,  the  net  thrust  is  also  the  equal  and 
opposite  of  the  net  change  of  momentum  of  the  internal  stream. 

Let  us  consider  now  three  streams  of  airflow:  one  over  the  upper 
surface  of  the  wing,  another  th&  internal  stream,  and  third,  the  flow 
over  the  lower  surface  of  the  wing.  The  internal  stream  after  passing 
through  an  oblique  shock  attached  to  the  leading  edge  of  the  upper  wing 
surface  and  through  a  normal  shock  just  upstream  of  the  leading  edge  of 
the  lower  wing  surface  arrives  at  the  entrance  of  the  subsonic  diffuser 
with  but  small  total  pressure  loss.  The  internal  flow  enters  the  subsonic 
diffuser  with  a  disturbance  on  the  upper  surface  due  to  the  boundary  layer 
and  the  normal  shock,  and  a  disturbance  on  the  lower  surface  due  to  the 
vorticity  of  a  curved  normal  shock.  The  process  of  subsonic  diffusion 
has  these  notable  features:  the  aforementioned  disturbances  at  the  inlet, 
the  low  Reynolds  number  of  the  flow  and  the  low  divergence  angle  and  long 
lengths  of  the:  diffuser  channel.  The  divergence  angle  and  lenth  of  the 
diffuser  channel  ia  the  result  of  the  compromise  in  the  design  for  maximum 
wing  loading  between  wing  profile  drag  and  diffuser  pressure  loss  on  the 
one  hand,  and  increased  lift  and  decreased  wing  pressure  drag  on  the  other 
hand.  The  wing  profile  drag,  the  diffuser  pressure  loss  and  the  wing  lift 
increase  with  diffuser  length;  the  wing  pressure  drag  decreases  as  the  diffuser 
section  divergence  angle  decreases.  The  internal  stream  arrives  at  the 
entrance  of  the  heat  transfer  section  with  vorticity  and  turbulence  due  to 
upstream  flow  conditions  and  including,  in  particular,  the  boundary  layer 
effects  in  the  diffuser.  It  is  expected  that  due  to  the:  low  Reynolds  number 
of  the  diffuser  flow  and  the  low  divergence  angle  of  the  diffuser  channel 
that  separation  will  not  occur. 
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In  the  heat  transfer  section  there  is  pressure  loss  due  to  heat 
transfer ,  somewhat  according  to  the  Reynolds  analogy.  There  is  also  a 
pressure  loss  due  to  heating,  accordingly  as  the  flow  is  accelerated  with 
decreasing  density.  Both  of  these  effects  increase  sharply  with 
increase  of  Mach  number.  The  Mach  number  of  the  heat  transfer  section 
flow  as  a  function  of  distance  along  the  heat  transfer  section  is 
critically  dependent  on  the  diffuser  section  outlet  Mach  number. 

The  diffuser  section  outlet  Mach  number  is  most  critically  dependent 
in  this  case,  on  the  diffuser  outlet  to  inlet  area  ratio.  This  ratio 
of  areas  also  effects  the  divergence  angle  and  length  of  the  diffuser 
section.  Besides  the  pressure  loss  in  the  heat  transfer  section  at  the 
heat  transfer  foils  due  to  heat  transfer  and  accelerated  flow,  there  is 
the  pressure  loss  at  the  channel  walls.  This  pressure  loss  also  increases 
with  the  local  Mach  number. 

The  rate  at  which  the  Mach  number  increases  along  the  heat  transfer 
channel  and  hence  the  rate  at  which  pressure  loss  occurs  with  friction 
and  heat  addition  depends  critically  on  the  length  of  the  heat  transfer 
foils  and  the  divergence  angle  of  the  heat  transfer  section.  The  effect 
of  the  divergence  angle  is  that  of  subsonic  diffusion.  The  shorter  the 
heat  transfer  foil,  the  higher  the  heat  transfer  coefficient  but  also 
the  higher  the  drag  coefficient.  There  exists,  however,  for  each  design 
case  an  optimum  combination  of  heat  transfer  channel  divergence  and  heat 
transfer  foil  length  for  the  compromise  of  pressure  drag  in  the  external 
stream  and  pressure  loss  in  the  internal  stream  with  heat  transfer  in  the 
internal  stream. 

The  optimum  length  of  the  heat  transfer  section  is  reached  when  with 
the  increase  of  the  internal  stream  Mach  number  and  with  the  increase  of 
the  internal  stream  and  heat  transfer-^oi-1  temperatures  the  rates  of 
internal  stream  pressure  loss  and  rate  of  radiation  loss  from  the  aircraft 
are  so  large  that  no  increase  of  wing  loading  occurs  with  increase  of  heat 
transfer  section  lengths. 

The  areas  of  the  jet  nozzle  throat  and  outlet  are  determined  by  the 
stream  conditions  at  the  entrance  of  the  jet  nozzle  throat  and  by  the 
requirement  that  the  static  pressure  in  the  jet  at  the  outlet  be  equal 
to  the  atmospheric  pressure. 
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This  requirement,  to  a  good  approximation  for  the  condition  of  the  solar 
powered  aircraft  jet,  determines  the  nozzle  areas  for  maximum  thrust. 

The  length  of  the  jet  nozzle  section  is  determined  for  optimum  wing 
loading  by  a  compromise  of  wing  pressure  and  profile  drags  with  the 
internal  stream  pressure  loss.  The  dominant  effect  is  sufficient  length 
and  hence  sufficiently  small  divergence  angle  that  the  pressure  drag  in 
the  external  stream  is  small.  It  appears  as  is  demonstrated  in  the  section 
on  the  results  of  design  and  performance  calculation  that  the  need  for  small 
pressure  drgg  in  the  external  stream  is  so  severe  that  the  divergence  angles 
of  both  the  subsonic  diffuser  and  the  jet  nozzle  sections  are  far  smaller 
than  optimum  for  the  internal  stream  flow. 

The  problem  of  the  design  and  performance  of  the  wing  in  the  external 
stream  for  the  upper  surface  is  familiar  and  likewise  for  the  lower  surface 
except  for  the  spillage  flow  from  the  subsonic  diffuser  at  the  leading  edge. 

A  detailed  account  is  given  in  the  following  of  the  derivation  with 
assumptions  and  approximations,  of  the  method  of  design  and  performance 
analysis. 
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5.  THE  PLAN  AND  BASIS  FOR  THE  AIRCRAFT  PERFORMANCE  CALCULATIONS 

The  aim  of  this  analysis  is  to  determine  a  valid  estimate  of  the 
performance  of  a  solar  powered  aircraft.  The  simplifying  assumptions 
and  approximation*  made  in  this  analysis  are  designed  for  either  negligible 
error  or  for  small  over- estimation  of  performance.  The  performance  est¬ 
imated  will  be  the  performance  at  best.  The  purpose  of  the  performance 
calculations  is  to  provide  a  basis  for  a  definitive  statement  of  the 
feasibility  of  a  solar  powered  aircraft. 

The  configuration  on  which  this  analysis  is  based  is  shown  in  Fig.  1 
and  is  discussed  in  Section  2.  The  flight  conditions  of  the  aircraft 
for  this  analysis  are  shown  on  Table  I  and  are  discussed  in  Section  3. 

The  processes  by  which  the  aircraft  flies  are  discussed  in  elementary 
terms  in  Section  4. 

For  convenience  the  aircraft  is  considered  in  three  major  parts: 
the  inlet  diffuser  from  Station  1  to  Station  2,  the  radiation  collector- 
heat  transfer  section  from  Station  2  to  Station  3 ,  and  finally,  the  jet 
nozzle  section  from  Station  3  to  Station  4.  The  diffuser  section  is 
considered  in  two  parts:  a  supersonic  diffuser  from  Station  1  to  Station  It 
and  a  subsonic  part  from  Station  It  to  Station  2. 

The  physical  basis  and  major  assumptions  of  the  aircraft  performance 
analysis  are  as  follows: 

1.  The  aircraft  basic  configuration 

2.  The  aircraft  flight  conditions 

3.  Estimated  pressure  recovery  ratios  for  the  supersonic  and 
subsonic  diffusers. 

4.  The  averages  of  the  flow  conditions  along  the  collector- transfer 
section  are  computed  as  if  the  flow  were  one -dimensional . 

5.  The  pressure  loss  and  heat  transfer  at  the  heat  transfer  foils 
are  estimated  by  flat  plate  laminar  boundary  theory. 
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6.  The  pressure  loss  and  heat  transfer  at  the  internal  channel 
walls  are  estimated  by  assuming  the  boundary  layers  to  be 
just  at  transition  from  laminar  to  turbulent. 

7.  The  profile  drag  and  beat  transfer  in  the  external  stream  at 
the  wing  upper  and  lower  surfaces  are  estimated  by  assuming 
the  boundary  layers  to  be  laminar „ 

8.  The  pressure  drag  of  the  wing  in  the  external  stream  is 
computed  by  the  Busemann  theory. 

9.  A  simplified  model  is  adopted  for  the  analysis  of  the  radiant 
heat  transfer  to  and  from  the  heat  transfer  foils.  By  means 

of  this  theory  and  the  optical  properties  of  the  various  surfaces 
involved  a  performance  method  is  evolved  for  the  estimation  of 
the  radiant  heat  transfer. 

10.  Estimated  pressure  recovery  ratio  for  the  jet  nozzle. 

These  physical  bases  and  major  assumptions  are  adopted  for  the  purpose 
of  estimating  the  best  performance  of  the  aircraft  which  may  be  expected. 

The  critical  item  of  performance  sought  is  the  attainable  wing  loading 

(^/xc)  a 

The  performance  analysis  method  is  derived  in  Section  11. 

The  performance  analysis  sequence  of  computations  is  given  in  Section  12. 
Brief ly,  the  analysis  procedure  is: 

1.  Choose  the  input  data 

Z  *  flight  altitude 

y^  =  inlet  diffuser  inlet  height 

b  =*  percent  of  span  to  be  used  as  diffuser  width 

R  =  diffuser  pressure  recovery  ratio 

pd 

R  »  jet  pressure  recovery  ratio 

“  collector-transfer  section  inlet  Mach  number 
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y/y2  “  collector- transfer  section  divergence  function 

xf/y  «  heat  transfer  foil  chord  length  ratio 


and  the  quantitites  E£r,  E  ,  E.,,  E  E.*,  E  which  refer 
to  the  radiation  transmitting,  reflecting,  absorbing  and  emitting 
properties  of  the  heat  transfer  foils  and  the  inner  and  outer  surfaces 
of  the  aircrafts 

2.  Proceed  with  the  operations  of  the  analysis  sequence  in  Section  12, 
These  operations  yield,  besides  much  detailed  information,  the  dimensions 
Jy  7^  and  xc  as  well  as  the  dimensional  ratios  i^/x^,  ^2^xc; 
all  with  the  performance  results:  the  excess  thrust  coefficient 


F 


-  C, 


and  the  attainable  wing  loading  (w/x 


The 


dimensions  and  dimensional  ratios  are  determined  for  the  maximum  value  of 
the  wing  loading 

3.  Repeat  step  (2)  for  various  values  of  y/y 2 ,  x^/y,  and  h 

until  the  set  of  these  is  determined  for  maxi urn  wing  loading  (w/x  1  a. 

4.  Repeat  steps  (2)  and  (3)  for  various  altitudes  until  the 
altitude  is  determined  for  the  maximum  wing  loading,  ^w/xj^ 

The  validity  of  these  results  of  performance  computation  depend  critically 
on  the  validity  of  the  data:  the  diffuser  and  jet  recovery  ratios,  R  ^  and 
Rp_j,  the  radiant  energy  transfer  parameters  Eff>  E^u,  Euf>  anc* 

E  and  on  the  validity  of  the  physical  basis  and  assumptions  (4), (5),  (6) 
and  (7)  which  determine  the  estimates  of  heat  and  momentum  transfers  to  the 
internal  stream  in  the  collector-transfer  section  and  to  the  external  stream 
on  the  upper  and  lower  wing  surfaces . 
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The  pressure  loss  and  heat  transfer  at  the  internal  channel 
walls  are  estimated  by  assuming  the  boundary  layers  to  be 
just  at  transition  from  laminar  to  turbulent. 

7.  The  profile  drag  and  heat  transfer  in  the  external  stream  at 
the  wing  upper  and  lower  surfaces  are  estimated  by  assuming 
the  boundary  layers  to  be  laminar. 

8.  The  pressure  drag  of  the  wing  in  the  external  stream  is 
computed  by  the  Busemann  theory. 

9.  A  simplified  model  is  adopted  for  the  analysis  of  the  radiant 
heat  transfer  to  and  from  the  heat  transfer  foils.  By  means 

of  this  theory  and  the  optical  properties  of  the  various  surfaces 
involved  a  performance  method  is  evolved  for  the  estimation  of 
the  radiant  heat  transfer. 

10.  Estimated  pressure  recovery  ratio  for  the  jet  nozzle. 

These  physical  bases  and  major  assumptions  are  adopted  for  the  purpose 
of  estimating  the  best  performance  of  the  aircraft  which  may  be  expected. 

The  critical  item  of  performance  sought  is  the  attainable  wing  loading 

(w\)  • 

The  performance  analysis  method  is  derived  in  Section  11. 

The  , performance  analysis  sequence  of  computations  is  given  in  Section  12. 
Briefly,  the  analysis  procedure  is: 

1,  Choose  the  input  data 

Z  =  flight  altitude 

y^  =  inlet  diffuser  inlet  height 

b  **  percent  of  span  to  be  used  as  diffuser  width 

R  ,  «  diffuser  pressure  recovery  ratio 

pd 

R  «  jet  pressure  recovery  ratio 

“  collector- transfer  section  inlet  Mach  number 
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y/y2  «  collector- transfer  section  divergence  function 

x^/y  *  heat  transfer  foil  chord  length  ratio 


and  the  quantitites  to.,  E  *  E*-*  E  E-jC,  E  which  refer 

't  ff'  uu;  uf  *  IV  su 

to  the  radiation  transmitting,  reflecting,  absorbing  and  emitting 
properties  of  the  heat  transfer  foils  and  the  inner  and  outer  surfaces 
of  the  aircraft. 

2.  Proceed  with  the  operations  of  the  analysis  sequence  in  Section  12. 
These  operations  yield,  besides  much  detailed  information,  the  dimensions 
y«>  y~,  y,  and  X  as  well  as  the  dimensional  ratios  i1/x  ,  i0/x  ,  i  /x  , 
all  with  the  performance  results:  the  excess  thrust  coefficient 


f—  -  o' 

b*c  dJ 


and  the  attainable  wing  loading 


H 


The 


dimensions  and  dimensional  ratios  are  determined  for  the  maximum  value  of 
the  wing  loading 


nsional  ratj 
(w/xj  a. 


3.  Repeat  step  (2)  for  various  values  of  y/y^,  xf/y>  a**d  h 

until  the  set  of  these  is  determined  for  maximum  wing  loading  |w/x  j  a. 

4.  Repeat  steps  (2)  and  (3)  for  various  altitudes  until  the 
altitude  is  determined  for  the  maximum  wing  loading,  ^w/x^ 

The  validity  of  these  results  of  performance  computation  depend  critically 
on  the  validity  of  the  data:  the  diffuser  and  jet  recovery  ratios,  R  ^  and 
R  j,  the  radiant  energy  transfer  parameters  Eff,  E^,  Euf,  E^f,  and 

Egu  and  on  the  validity  of  the  physical  basis  and  assumptions  (4),  (5),  (6) 
and  (7)  which  determine  the  estimates  of  heat  and  momentum  transfers  to  the 
internal  stream  in  the  collector- transfer  section  and  to  the  external  stream 
on  the  upper  and  lower  wing  surfaces. 
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6.  PRIOR  DESIGN  AND  PERFORMANCE  ANALYSIS  OF  THE  SOLAR  POWERED  AIRCRAFT 
In  WADC  TN59-184,  the  final  report  of  the  first  phase  of  this  work, 
the  results  of  a  simpler  but  less  valid  analysis  are  presented.  These 
results  are  shown  in  Table  2.6-1  of  WADC  TN59-184  and  in  Fig „  2 
of  this  report.  The  design  and  performance  results  of  Table  2.G-1  of 
WADC  TN59-1S9  led  to  the  choice  of  the  parameters; 

M1  * 

M2  »  0.3 

%  “  °-6 


Pd 


0.03 

0.93 


for  the  design  of  the  optimum  solar  powered  aircraft;  the  results  of  the 
performance  analysis  of  designs  based  on  these  parameters  are  shown  in  Fig.  2 
The  performance  results  in  Fig. 2  are  shown  here  not  in  direct  support  of  a 
statement  of  the  feasibility  of  a  solar  powered  aircraft  but  as  an  indication 
of  the  optimum  altitude  for  a  very  large  solar  powered  aircraft.  The  choice 
of  the  altitude,  200,000  ft.  for  the  design  and  performance  analyses  of  this 
report  is  based  on  Fig.  2.  It  is  believed  that  the  trends  of  chord  length, 
altitude  and  wing  loading  in  Fig.  2  are  significant  but  that  the  magnitude 
of  the  attainable  wing  loadings  are  greatly  overestimated.  The  particular 
assumptions  of  the  prior  method  of  design  and  performance  analysis  which  led 
to  the  simple  computation  procedure  with  its  extensive  but  insufficiently  valid 
results  are: 

1.  The  heat  transfer  airfoil  temperature,  T^,  is  constant  along  the 
col lector- transfer  section  and  equal  to  1133°K. 

2.  The  air  temperature  at  the  outlet  of  the  collector  transfer 
section,  T^,  is  equal  to  Tf. 

3.  The  pressure  recovery  ratio  across  the  heat  transfer  section, 

P  /?  —  P.  is  estimated  as  if  the  Mach  n usher  were  constant 

s3  s2  ph 

along  the  heat  transfer  section. 

The  design  and  performance  method  of  this  report  was  devised  as  a  major 
step  toward  sufficient  validity. 
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7 .  RESULTS  OF  SOLAR  POWERED  AIRCRAFT  DESIGN  AND  PERFORMANCE  ANALYSIS 

The  input  data  are  shown  in  Table  II;  certain  design  and  performance 
results  are  listed  in  Table  III.  The  corresponding  aircraft  configurations 
are  indicated  in  Figs.  3a  and  3b.  A  brief  description  of  the  computational 
procedure  is  given  in  Section  5;  the  details  procedure  is  set  forth  in 
Section  12. 

The  optimum  aircraft  for  200,000  feet  altitude  and  diffuser  inlet 

height  of  20  feet  is  indicated  by  the  third  configuration  of  Case  XI;  the 

2 

attainable  wing  loading  is  .042  lbs/ft  for  a  chord  length  of  1700  feet. 

In  the  design  and  performance  study  of  this  report  indicated  by  Table  II, 
Table  III,  and  Figs.  3a  and  3b,  the  flight  altitude  and  the  diffuser  inlet 

the  Mach  Number  at  the 


V 


height  were  held  constant  and  the  quantities: 

diffuser  outlet,  d/dx  ( ’y/y the  divergence  of  the  collector- transfer  section; 

x^/y,  the  heat  transfer  foil  length;  x^,  the  wing  chord  length,  were  varied 

until  the  optimum  set  of  them  was  found.  It  is  interesting  to  note  from 

Tables  II  and  III  the  critical  dependence  of  the  wing  loading  w/x  upon 

c 

the  set  of  quantities:  d/dx  (y/y2)>  (x^/y)  and  xc *  There  is  needed 

the  repetition  of  this  procedure  for  the  same  diffuser  inlet  height  for 
other  altitudes  in  order  to  find  the  optimum  solar  powered  aircraft  according 
to  the  methods  of  this  report.  The  configurations  shown  in  Figs.  3a  and  3b 
are  surprisingly  slender.  It  is-  a  clear  result  of  the  prior  design  and 
performance  analyses  shown  in  Table  2.6-1  of  WADC  TN59-184  that  the  wing 
thickness  to  chord  length  ratio  for  optimum  configurations  varies  approximately 
inversely  at  the  atmospheric  density.  It  is  expected,  therefore,  that  the 
optimum  solar  powered  aircraft  according  to  the  methods  of  this  report  will, 
when  found,  have  higher  wing  thickness  ratio  than  shown  in  Figs.  3a  and  3b 
for  200,000  feet  flight  altitude. 

In  Figs.  4,5,6,7,8,9,10,11,12,13,14,15,16  for  Cases  I, II, III, IV, V, VI, 

VII, VIII, X, XI, XII, XIII, XIV  there  is  plotted  the  thrust,  F,  in  lbs.,  the 
Mach  number  M,  the  temperature  ratios: 

!* 

T* 

u 

T! 

T* 

s 

each  ratio  relative  to  Tg2>  the  internal  stream  stagnation  temperature  at  the 
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inlet  of  the  collector  transfer  section  and  th e  fractional  part  of  the 

incident  solar  energy  transferred  to  the  internal  stream,  all  as  functions 
of  the  distance  x  from  the  wing  leading  edge. 

A  survey  of  the  curves  in  Figs*  4,5,6,7,8,9,10,11,12,13,14,15,  and  16 

reveal  the  critical  dependence  of  F,  the  total  thrust,  upon  T*,  p*,  and 

M0  and  the  critical  relation  between  p*  and  .  Further  consideration  of 
*3  s  3 

these  curves  reveal  that  the  most  favorable  condition  for  maximum  (w/x^)a 
occurs  with  small  and  substantially  constant  with  increasing  x. 

This  condition  is  obtained  with  the  proper  set  of  quantities  d/dx  (y/y^)  and 
s^/y.  These  effects  are  due  to  the  extremely  critical  dependence  of  the 
total  thrust  on  the  process  of  convective  heat  transfer  and  pressure  loss  due 
to  friction  occuring  in  the  internal  stream  at  the  collector  transfer  foils 
and  the  extremely  critical  dependence  of  this  heat  transfer- pressure  loss 
process  on  the  parameters  d/dx  (y/y^)  and  (x^/y).  The  parameter  d/dx  (y/y^) 
is  also  critical  with  respect  to  drag  of  the  external  stream. 

In  Figs.  17,  18,  and  19  for  cases  VIII,  XI,  and  XII,  are  plotted  the 
terms  of  Eq.  (6)  for  the  energy  flux  to  the  internal  stream.  The  values 
of  0.915  indicates  that  91.5  percent  of  the  insolation  on  the  collector- 
transfer  section  is  absorbed  by  the  heat  transfer  foils;  it  is  estimated 
that  loss  of  foil  emission  by  transmission  through  the  upper  wing  surface  is 
negligible.  The  loss  of  energy  at  wing  upper  and  lower  surfaces  is  principally 
by  emission;  the  loss  there  by  convection  to  the  external  stream  is  small. 

The  heat  transfer  to  the  wing  surfaces  is  principally  by  radiative  exchange 
between  the  collector  transfer  foils  and  the  wing  surfaces;  the  heat  transfer 
by  convection  from  the  internal  stream  to  the  wing  surfaces  is  negligible. 

And,  of  course,  the  heat  transfer  from  the  collector- transfer  foils  to  the 
internal  stream  is  by  convection.  The  very  significant  effect  to  be  noted 
here,  is  that  a  major  loss  of  energy  occurs  by  radiation  from  the  collector- 
transfer  foils  to  the  wing  upper  and  lower  surface  and  hence  by  emission  to 
space . 
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8.  VALIDITY  OF  THE  RESULTS  OF  PERFORMANCE  ANALYSIS 

The  results  of  the  performance  analysis  are  summarized  in  Table  III. 

These  results  were  reached  by  wall  considered  methods  of  analysis  using 

such  experimental  results  as  were  available  and  applicable.  There  is 

much  to  be  desired,  however,  in  the  correspondence  of  both  the  simple 

models  used  as  a  basis  of  the  performance  analysis  and  the  elementary 

experimental  results  with  the  configurations  and  processes  of  a  real 

solar  powered  aircraft.  In  spite  of  these  doubts  there  remains  the 

results  of  Table  III.  The  most  significant  result  in  Table  III 

is  the  third  configuration  of  Case  XI;  the  attainable  wing  loading 

2 

estimated  is  0.042  Ibs/ft  for  a  chord  length  of  1700  ft.  flying  at  air 
altitude  of  200,000  ft.  It  is  suggested  that  better  performance  with 
shorter  chord  length  will  be  obtained  at  a  higher  altitude  less  than 
220,000  ft. 

Without  consideration  of  the  value  of  the  solar  powered  aircraft, 
the  results  shown  in  Table  III  are  judged  from  a  technical  point  of  view 
to  be  of  sufficient  validity  to  justify  the  following  steps: 

1.  More  Complete  Design  and  Performance  Analysis: 

a.  Revise  and  develop  in  detail  the  performance  methods 
of  this  report . 

b.  Seek  extensively  and  assiduously  experimental  results 
in  the  literature  for  the  support  of  the  input  data  of 
design  and  performance  analyses, 

c.  Accomplish  design  and  performance  analyses  over  a  sufficient 
range  of  the  parameters: 

Z;  the  flight  altitude;  y^,  the  diffuser  inlet 

height;  M^,  the  diffuser  outlet  Mach  number  d/dx,  (y/y^) 

the  divergence  of  the  collector  transfer  section. 

*  /y  ,  the  heat  transfer  foil  length  ratio 
x  ,  the  wing  chord  length. 


9.  THE  UNSYMMETRICAL  TWO-SHOCK  INLET  DIFFUSER 

The  unsymmetrical  two-shock  inlet  diffuser  is  described  and 
discussed  in  Refs.  2  and  3.  In  particular,  attention  is  directed  in 
Ref.  2  to  Section  7  and  Figs.  7.01  and  7.19.  A  sketch  of  the  unsym- 
metrical  two-shock  inlet  diffuser  with  its  shock  wave  configuration 
for  the  design  Mach  number  of  1.5  is  shown  in  Fig.  20.  The  stagnation 
pressure  ratio  across  the  two-shock  system  from  free  stream  M  =  1.50 
to  inlet  M*  0.860  is  0.9867.  The  stagnation  pressure  ratio  just  down¬ 
stream  of  the  normal  shock  is  slightly  less  than  0.9867 ,  due  to  the 
boundary  layer  on  the  upper  surface.  The  advantages  of  this  design  are; 

1.  There  is  no  swallowing  problem  since  the  entrance  section 
and  throat  are  identical. 

2.  In  the  supersonic  part  of  the  diffusion,  from  free  stream  to 
within  the  normal  shock,  there  is  boundary  layer  effect  on 
only  one  surface . 

3.  The  unsymmetrical  condition  due  to  the  boundary  layer  at  the 
normal  shock  on  the  upper  surface  insures  that  the  unsym¬ 
metrical  boundary  layer  thickening,  turbulence  and  separation 
effects,  which  occur  even  in  symmetrical  subsonic  diffusers, 
will  occur  in  this  case  on  the  upper  surface,  and  a  more 
favorable  boundary  layer  may  be  maintained  on  the  lower  surface. 

The  prospect  of  localizing  the  turbulence  and  separation  in  the 
subsonic  diffuser  on  the  upper  surface  suggests  starting  a  diagonal 
array  of  heat  transfer  foils  from  the  upper  surface  and  from  within 
the  subsonic  diffuser.  Thus  the  turbulence  occurring  in  this  region 
will  be  utilized  for  heat  transfer  and  the  spread  of  large  scale  eddies 
will  be  reduced  by  means  of  the  heat  transfer  foils  acting  as  vanes. 
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10.  ELEMENTARY  CONSIDERATIONS  OF  TOTAL  THRUST  AND  REYNOLDS  ANALOGY 


10.1  Discussion  of  tha  Total  Thrust 

The  Equation  124  estimates  the  total  thrust  per  unit  span 
for  the  solar  powered  aircraft.  Let  us  derive  an  expression  for  the  total 
thrust  per  unit  span  as  in  Section  11.3  but  in  this  case: 

1.  Utilizing  only  a  fraction  b  of  the  span  for  diffuser 

inlet,  and 

2.  Estimating  the  temperature  R,^  as 


1  + 


t8S  %  xc 
pl^lbylcpT 


1  + 


g  n  x 
s  *s  c 


(Nn).  b(N  ).  RT.. 
v  R  i  v  p^l  1 


a) 


where  n  is  the  net  fraction  of  the  incident  insolation  which  is  delivered 

's 

to  the  internal  stream  and  is  not  lost.  Then, 


F  *  2gby 


VC 


i  + 


gs  %  xc 


V]  -  (-1 


where 


o  <  b  4  1 


b<wyiRTi  j 


(Y-l) 


*p/  Y 


-  1 


(2) 


Thus,  F  ia  a  function  of  the  design  parameters  b,  x^ ,  t|g;  it  is  a 
monotonically  increasing  function  with  all  these  parameters.  The  use  of 
b  ■  1  rather  than  a  lower  value  will  determine  that  the  internal  stream 
temperature  will  be  lower.  Also,  due  to  the  lower  stream  temperature  the 
foil  temperature  will  be  lower  and  hence  tj  will  be  higher  since  the  loss 
of  the  heat  delivered  to  the  foils  is  chiefly  by  radiation  to  the  upper  and 
lower  wing  surfaces.  Therefore,  except  as  R^,  the  pressure  recovery  ratio  is 
related  to  b,  and  x^,  the  total  thrust  per  unit  span  is  maximum  for  b  *1. 

Aa  we  shall  see  in  Section  10.2,  the  pressure  recovery  ratio  is,  however,  related  to 
b>  *\s  c 
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10.2  Relation  Between  Pressure  Loss  and  Heat  Transfer 


According  to  Eq.  63  of  Ref.  7  in  the  notation  of  this  report, 


Of 

R(T£  -  T  ) 
t  e 


cd*  Re  r 

f  xf 


(3) 


This  equation  expresses  the  Reynold’s  analogy  relation  between  skin  friction 
and  heat  transfer  for  a  flat  plate  with  constant  surface  temperature..  The 
pressure  loss  per  unit  length  is 


pu2cdf 

y 


(4) 


Then,  from  (3)  and  (4)  with 
P  =  pRT 


*«,  -  V  Vp  1/3 


(5) 


The  ratio  of  the  rate  of  stagnation  pressure  loss  to  the  rate  of  heat 
transfer  in  the  collector- transfer  duct  increases  as  the  square  of  the 
local  Mach  number,,  inversely  as  the  temperature  difference  of  the  heat 
transfer  foils  and  the  internal  stream  and  inversely  as  the  local  Reynolds 
number  referred  to  stream  height.  Thus,  reduced  airflow  by  means  of  b  <  1 
as  discussed  in  Section  10.1  may  allow  such  reduction  in  stagnation  pressure 
loss  as  to  yield  increased  thrust  and  finally  increased  wing  loading.  This 
matter  is  judged  to  be  so  involved,  however,  as  to  require  the  inclusion 
of  b  with  the  parameters  to  be  varied  in  seeking  the  oDtimum  design.  Since, 
M2  1 

—  in  (5)  varies  approximately  as  it  is  clear  that  the  maximum  cross- 

ey  y 

section  of  the  internal  stream  is  consistent  with  the  external  drag  is 
desirable  for  the  reduction  of  the  pressure  loss  in  the  internal  stream. 
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11. 


DEVELOPMENT  OF  METHOD  OF  ANALYSIS 

11-1  Performance  of  the  Radiation  Collector  Heat  Transfer  Section 


11.1.1  The  Equations  of  Energy  Flux  in  the  Collector -Tr aha fer 
Section 


where 


H  = 
*0  S 


C 


int 


The  energy  flux  to  the  internal  stream  is 


42 

dx 


n  q, 

lo  s 


-  R  -  C 


ext 


+  C 


int 


(6) 


energy  flux  from  the  sun 

energy  flux  from  the  sun  absorbed  by  the  collector- 

transfer  section 

energy  flux  lost  by  radiation 

energy  flux  lost  by  convection  to  the  external  stream 
energy  flux  by  convection  from  the  collector-transfer 
foils  to  the  internal  stream 

energy  flux  by  convection  from  the  col lector -transfer 
wall  surfaces  to  the  internal  stream. 


and  where  the  energy  fluxes  rj  <1  ,  R,  C  ,  C,-,  and  C,  are  defined 

o  s  ext  t  int 

as  energy  per  unit  time  per  unit  area  normal  to  the  solar  beam.  These 
energy  fluxes  are  functions  of  position  along  the  wing  chord. 

The  condition  that  the  net  flux  of  energy  to  the 
upper  and  lower  wing  surfaces  are  each  zero  yields  the  equations: 


^-(f-m)  (u-*f)  Cint  u^  ^ext  +  u->  +  u— ^ 

f  7  u  '  u  S  u  u  f 

l)"\1(l->f)"Cint1  =  C  +  %1  +  V-  (8) 

where 


q\f(f-*u) 


Collector-transfer  surface  thermal  radiation  energy 
absorbed  by  the  wing  upper  surface. 
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(u-»f)  = 
u 


(1— >f  )  = 


f> 


<W 

s 


X:U-» 


V- 


Wing  upper  surface  thermal  radiation  energy  absorbed 
by  the  collector-transfer  surface. 

Collector -transfer  surface  thermal  radiation  energy 
absorbed  by  the  wing  lower  surface. 

Wing  lower  surface  thermal  radiation  energy  absorbed 
by  the  collector-transfer  surface. 

Solar  radiation  absorbed  by  the  collector-transfer 
surface . 

Solar  radiation  absorbed  by  the  wing  upper  surface. 

Solar  radiation  reflected  externally  at  the  upper  wing 
surface . 

Wing  upper  surface  thermal  radiation  emitted  externally. 

Wing  lower  surface  thermal  radiation  emitted  externally. 

Collector- transfer  surface  thermal  radiation  energy 
transmitted  by  the  wing  upper  surface. 

Collector- transfer  surface  thermal  radiation  energy 
transmitted  by  the  wing  lower  surface. 


11.1.2  Radiant  Energy  Transfer  in  the  Collector-Transfer 
Section 

The  arrangement  of  radiation  collector-heat  transfer 


foils  under  consideration  for  the  solar  powered  aircraft  is  as  shown  in 


Fig.  21.  The  analysis  of  the  radiant  energy  transfer  in  this  configura¬ 
tion  has  proved  very  difficult.  Therefore,  in  order  to  obtain  a  rational 
and  an  "at  best"  estimate  of  the  radiant  energy  transfer,  assume  the 
collector- transfer  foils  to  be  equivalent  to  a  plane  sheet  between  the 


upper  and  lower  surfaces  of  the  wing  as  in  Fig.  22,  where 


«,  p,  OC,T  =  emissivity,  reflectance,  absorptance  and  transmittance, 
respectively. 

s,  u,  f,l  *  subscripts  referring  to  the  sun,  wing  upper  surface, 
collector-transfer  surface  and  wing  lower  surface, 
respectively. 
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XsAuAf  Ai  "  subscripts  referring  to  solar  radiation,  thermal 

radiation  of  wing  upper  surface,  collector- transfer 
surface  and  wing  lower  surface,  respectively. 


In  general, 


Xij 


(9) 


V  -  %)  +  v  -  1  <l0> 

Consider  first,  the  radiant  energy  exchange  between  two 
parallel  infinite  plane  surfaces  as  in  Fig.  23.  According  to  Ref.  4, 
Article  14.2,  p.  403,  the  heat  flux  lost  by  surface  1  as  radiation  is 
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and  the  radiation  lost  by  surface  2  as  radiation  is 

1  l 

\2(2  -  1)  1  P\22P\21  X2  2 


(12) 


Consider  now,  the  "green  house"  effect  between  two 
ite  plane  s 

retained  by  surface  2  is 


parallel  infinite  plane  surfaces  as  in  Fig.  24.  The  energy  flux 
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Now,  for  the  collector- transfer  section,  by  use  of  Eqs.  11  and  12 
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and  by  use  of  Eqs.  13  and  14 
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(17) 
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Assume 


q.  ■»  0  *  q.  - 
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(24) 


By  definition  of  tJoqg  with  Eqs.  19  and  20 

qXf  +qxu 

s  s 


s 

Tx  u  <L  '  px  f> 

s _ s 

1  ■  px  U  Px  f 

S  8 


(25) 


+  °X  u 

S 


1  + 


Xu  H Kf 

s _ s 

1  -  px  U  Px  f 

s  s 


and  by  definition  of  R  with  Eqs.  22,  23  and  24 
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11.1.3  The  Temperature  Change  Along  the  Internal  Stream 

It  is  assumed  that  the  net  energy  flow  to  the  internal 
airstream  occurs  uniformly  across  the  stream.  Hence, 
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where 

y«p  ■  y2u2P2  (28) 

11.1.4  One  Dimensional  Channel  Flow  in  the  Radiation 

Collector-Heat  Transfer  Section 

The  approximation  of  one  dimensional  channel  flow 
is  better  as  the  arrangement  of  the  collector- transfer  foils  tends  toward 
a  uniform  density  of  foils  across  each  cross-section  thus  tending  toward 
uniform  heat  flux  to  the  stream  and  uniform  flow  conditions  across  each 
section. 

For  one  dimensional  channel  flow  with  heat  addition, 
friction  and  area  change,  according  to  Ref.  5: 


^  VM2  r  dA 

j  ttw  1^7 


(29) 


(30) 


where  P  ,  T  t  A  ,  M0  are  stream  conditions  at  Station  2.  The  subscript  s 
S2  s2  1  1 

refers  to  stagnation  conditions. 


(31) 

(32) 


where  P  and  T  are  the  stream  conditions  of  pressure  and  temperature. 
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M  3  Mach  number 

2  x  Area  of  stream  cross-section 

m  sa  - - * - 

Wetted  Perimeter  of  stream  cross-section 


m  *  hydraulic  radius 

f  «  channel  flow  friction  factor,  such  that 

(s)  -  zt  ■ 

\dx/  due  to  friction 


(33) 


(34) 


Equation  36  holds  with  or  without  either  or  both  friction 
and  heat  addition.,  When  A  and  f/2m  are  known  functions  of  x,  the 
equations  36  and  37  may  be  solved  for  Pg  and  M  as  functions  of  x. 

Equation  29  may  be  written  in  a  form  more  convenient  for 


use, as 
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where 


where  A/A*  is  tabulated  in  Ref.  6. 
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(36) 


11.1.5  Heat  Flux  by  Convection  from  the  Collector-Transfer 
Foils  to  the  Internal  Stream 

Assume  the  heat  transfer  from  the  foils  to  the  airstream 
is  for  each  foil  as  if  it  were  in  conventional  flat  plate  laminar  flow.  Then 
according  to  Eq.  53  of  Ref.  7  and  Eq.  14.68  of  Ref.  8: 
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where 
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*  Net  heat  flow  from  a  collector-transfer  foil  to  internal 
stream  per  unit  of  wing  span 
»  Chord  length  of  the  collector- transfer  foils 
«  Ratio  of  specific  heats 

ji,  T,  M  ■  Thermal  conductivity ,  coefficient  of  viscosity,  temperature 
and  Mach  number  of  the  internal  stream,  all  assumed  constant  across 
each  stream  cross-section. 


11.1.6  Momentum  Transfer  between  Collector-Transfer  Foils 
and  the  Internal  Stream 

Assuming  conventional  flat  plate  laminar  flow  on  the 
collector- transfer  foils,  then  according  to  Eq.  58a  of  Ref.  8: 


For  two  sides  of  one  foil 
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Hence 
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11.1.7  Estimation  of  the  Momentum  and  Convective  Heat  Transfer 

of  the  Collector-Transfer  Section  Walls  to  the  Internal  Stream 
On  the  assumption  that  the  drag  and  heat  addition  are  uni¬ 
formly  distributed  over  the  stream  cross-sections  and  on  the  fact  that  the 
momentum  transfer  from  the  walls  is  small  in  comparison  with  that  from  the  collector- 
transfer  foils,  the  velocity  distribution  across  the  stream  will  tend  to  be  uniform 
except  near  the  walls.  Therefore,  it  is  assumed  that  the  boundary  layers  are  just 
fully  developed  turbulent  throughout  the  channel  and  corresponding  to  the  point, 

N  =  —  m  4000  (46) 
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100  X  -  4 

on  Fig,  20.17,  page  418  of  Ref.  8.  Then 


X^2f  .04  (47) 

And  further,  according  to  empirical  formulas  of  Dittus  and  Boelter 
cited  on  page  211  of  Ref.  4 
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h  *  film  heat  transfer  coefficient  per  unit  area  of 
channel  surface.  Hence,  for 
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2m  m 


16k 

2m 

4  x  Area 


Wetted  Perimeter 
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w  »  width  of  duct 

The  heat  flux  by  convection  from  the  upper  and  lower  surfaces 
of  the  wing  to  the  internal  stream  are 
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where 
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temperature  of  the  upper  and  lower  wing  surfaces , 
respectively 
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1 1 . 1 . 8  Heat  Transfer  from  Upper  and  Lower  Wing 

Surfaces  to  the  External  Stream 

According  to  Eq.  50  on  page  557  of  Ref.  7 
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where 

q  »  local  rate  of  heat  transfer,  per  unit  area 
subscripts 

oo  *=■  condition  outside  the  boundary  layer 
w  ®  condition  at  the  wing  surface 

e  **  equilibrium  conditions  for  all  insulated  surface 
superscript 

*  *  dimensionless  quantity 
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Tw  *  wing  surface  average  temperature 
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where  t  (**)  may  be  approximated^  by  a  polynomial 

n 

t(x*)  -  ^  an  x*n 

a  «  t^  (Q)/n! 


J 


(59) 


T  -  T 

-w_e  (60) 


(61) 


where 


0 

1 

2 


*>> 

-0,5915 

-0,9775 

-1.1949 


Assume  skin  temperature  distribution  to  be  linear  with 
distance  from  the  leading  edge,  where  it  has  the  free  stream  recovery 
temperature 
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to  a  maximum  at  the  wing  trailing  edge.  Then, 
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Hence,  in  the  notation  of  this  report 
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11.1.9  The  Performance  Equations  of  the  Collector- 
Transfer  Section 

Collection  of  Results 
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By  Sutherlandts  formula 
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according 

From 


to  Ref.  7 

Eq.  6  with  Eqs.  27,  28,  26,  65  and  66,  37,  52  and  53 
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From  Eq.  8  with  Eqs.  17,  18,  23,  24,  52  and  65 


From  Section  11.1.4,  Eq.  32  and  Section  11. 1.5,  Eq.  39 


In  the  last  term  of  the  first  part  of  Eq.  95  put 
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The  approximations,  Eqs.  98,  99,  and  100,  achieve  an  important 
simplification  by  rendering  Eqs,  103,  104  and  10 5  independent  of 
Mach  number. 

From  Eqs*  29,  30 
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and  from  Eq*.  44,  47,  and  94 
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The  solution  of  Eqs.  103,  104  and  105  with 

Eq.  51  determine  T  ,  T  ,  T-  ,  and  Tr  as  functions  of  x.  Then 
s 1  u*  1 9  f 

Eq.  30  determines  M  as  a  function  of  x;  finally,  Eq.  29  determined 
Pg  also  as  a  function  of  x.  These  results  complete  the  analysis  of  the 
collector- transfer  section  except  with  respect  to  the  drag  of  the  external 
stream  and  with  respect  to  the  optimum  length  of  the  collector-transfer 
section. 

11.2  Performance  of  the  Inlet  Diffuser 

The  inlet  diffuser  configuration  considered  is  shown  in  Fig.  20 
and  is  discussed  in  Ref,  9.  The  critical  performance  parameter  of  the 
diffuser  is  the  stagnation  pressure  ratio  between  the  stations  1  and  2, 
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In  this  analysis,  is  assumed  to  be  a  known  constant.  The  flow  is 

assumed  adiabatic  from  station  1  to  station  2,  thus 
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By  means  of  Eq.  35 
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The  geometry  of  the  asymmetrical  inlet  and  the  subsonic 
diffuser  yield 
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11.3  Design  of  the  Jet  Nozzle  and  the  Determination  of  Net  Thrust 
For  the  flow  from  station  3  to  the  jet  outlet  at  station  4 
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The  net  thrust  per  unit  span  due  to  the  change  of  momentum 
of  the  internal  stream  from  station  1  to  station  4  is 
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For  convenience  and  near  optimum  performance 

P4  -  Pj  (123) 

Hence,  by  means  of  Eqs.  114,  116,  121,  122,  and  123,  Eq.  120  becomes 
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11,4  Estimation  of  the  Drag  of  the  Wing  External  Surface 

The  drag  of  the  wing  external  surface  is  estimated  as  the 
sum  of  a  pressure  drag  and  a  boundary  layer  friction.  The  pressure 
drag  coefficient  according  to  the  Busemann  theory  as  given  in  Ref.  10, 
Chapter  7 ,  Eqs*  266  and  267  is 
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and  where  according  to  Eqs.  255  and  260  of  Ref.  10  for  the  outside  surfaces 
of  the  particular  case  of  a  symmetrical  profile  shown  in  Fig.  1 
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(134) 


Assuming  laminar  boundary  layers,  as  on  a  flat  plate,  on  both 
the  upper  and  lower  surfaces  of  the  wing. 
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The  total  wing  drag  per  unit  span  is 
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11.5  Determination  of  Lengths  of  Inlet  Diffuser,  Collector- 

Transfer  and  Jet  Nozzle  Section  for  Maximum  w/x 
- : - 1 — c 

In  a  condition  of  steady  flight 

F  *  qxc  CD  (138) 

W  =  qXc  CL  (139) 

Where  F  and  w  are  the  thrust  available  and  the  aircraft 
weight  both  per  unit  span  and  where,  from  Eq.  137 
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On  eliminating  C  between  Eqs.  138,  139  and  140  we  have 
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We  seek  now  the  set  of  1^,  1^,  1^  which  maximize  the  quantity 
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to  Eq.  141  and  the  restriction,  Eq.  143.  For  this  purpose  put 
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and  from  Eq.  135 


From  Eqs.  145,  147,  149  and  151 
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12.  PERFORMANCE  ANALYSIS  SEQUENCE 

12.1  Throughout  the  analysis  the  following  quantities  are  assumed 
constant: 

qg  =  $2.4  ft  lb/ft^  sec 

y  *  1.40 

C  ,  air  «  6.01  x  10^  ft  lb/slug  °R 

P  -12  ,  2  o  4 

a  =  374  x  10  ft  lb/ft  sec  F 


12.2  The  properties  of*  air 

2 

|i  =  viscosity  (lb  sec/ft  ) 

k  “  thermal  conducting  (ft  lb/ft  sec  °R) 

N  *  Prandtl  number 

P 

are  plotted  from  Ref.  11  in  Figs.  25,  26,  and  27. 

12.3  Choose  the  following  sets  of  quantities; 


a.  Z,  Mp  y1 
b-  M2’  RPd'  W/,yl 

c.  Y/ y^,  xf/y  both  as  functions  of  distance  along  the 
collector- transfer  section 

d.  x  /x  ,  r|  ,  E’,  E  ,  E- 1  ,  E  ,  E  ,  E 

p  V  lo7  ff7  uu;  ll7  uf7  If7  su 

12.4  From  the  A  RPC  Standard  Atmosphere  obtain  a^,  jx  ,,  p^,  T^ 
and  corresponding  to  the  altitude  Z. 


12.5  Compute  u^,  (NR/x)^,  Tg 
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12.6  Compute  y^,  T2,  u2,  P^Pp  P2 
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where  m  (M)  =  1/(A/A*)  and  (A/A*)  is  given  in  Table  I  of  Ref.  6 
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12.7 

air  in  Figs. 


Obtain  \i^,  and  (Np)2 

25,  26 ,  and  27. 
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The  data  required  for  the  solution  of  the  equations  of 


Article  9  are: 


the  constants  t)  ,  216/T  ,  C  9  Cn,  Cy  •  Cc, 

xo7  s’  ^7  2 7  3 7  k7  5} 


c6,  C7-  Cg,  C9;  C10,  xd,  (w/yi),  yjy2,  the  plots  of  k,  Np  and  the 
functions 
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In  order  to  integrate  the  first  equation  in  186  ,  it  is  necessary 

*  *  * 

first  to  determine  a  solution  Tu  ,  and  for  the  last  equation 
of  186  with  187  and  188  .  For  the  latter  purpose,  for  given  Tg 
at  x  and  a  trial  value  of  T^*,  compute 
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Then,  with  the  available  data  compute 
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and  iterate  until  AT^  =  0 


Now,  compute  with  the  available  data 
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and  iterate  until  AT  =0.  Finally,  iterate  the  process  starting 
u  * 

with  190  and  ending  with  200  until  T  corrected  obtained  from 

*  “ 

198  is  the  same  as  the  Ty  introduced  as  a  trial  value  in  190  . 
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In  case  the  procedure  involving  the  set  of  equations 
192  -  200  does  not  converge,  then  the  following  procedure 

may  be  more  convenient.  Instead  of  194  and  199  use  iteratively 
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dTf* 
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*3 


These  processes  yield  the  set 


*  *  *  * 
x>  Tx  ,  Tu,  rv  Tf,  f2 


which  for  given  Tg  at  x  satisfy  the  set  of  equations:  the  last  o 
187  and  188. 

Plots  of  the  function 


1 


T  + 


216 


,1/4 


is  a  valuable  aid  in  the  computation  according  to  Article  12.9. 


12.10  Compute 
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=  T  '  \ 

!(i  +  1)  S(i)  J 
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i  4-  1 


dx 


"10 


by  means  of  "predictors”  and  "correctors".  The  "predictors"  and 
appropriate  to  the  points 


X1  "  Xd'  x2  “  X 


1  4-  §x,  x3  =  x9 


+  6x 


x~  +  c)x,  x.  .  _  *  x,  4*  §* 
3  x  4-  I  i 


(205) 

(206) 
f  186  , 

(207) 

(208) 

"correctors" 

(209) 
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are  listed  below.  The  value  of  T  >.  _  .  *  is  first  estimated  by  means  of 

x(i  -r  1}  y 

the  "predictors.”  Then,  with  this  as  the  given  value  of  T^  in  Article  12.9, 
the  set  201,  is  computed  and  in  particular  +  ^  is  computed  which 

substitued  in  the  "corrector"  yields  a  corrected  value  of  Tg .  Iteration  of 
this  process  will  yield  a  set  201  which  satisfy  208,  209,  and  210. 

At  the  first  point 


x  »  x  ,  T 
1  d7  s 
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u  ^first  trial 


lf  *Tf  ^  first  trial  "  1 


(210) 


At  the  second  point 


X2  “  X1  +  ^X 


^Tu  )  first  trial  *  ^u^  7  ^f  ^  first  trial  ^ 
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'Correctors 
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and  where  T  and  Tr  are  known  functions  of  x  tabulated  as  results 
s  f 

of  Articles  12.9  and  12.10 

The  input  data  for  Article  12.11  are  besides  Te  and  Tf  ,  the 


constants 


(221) 


and  the  functions 


(222) 


Denote  the  quantity  of  the  left  of  218  by  I  and  denote  the 
integrand  on  the  right  by  V  .  At  the  first  point  x^  =  x^ 


M  = 


1  = 


(223) 


7t’he  method  of  solution  proceeds  by  estimating  at  x  =  x^  a  trial 

value  of  I  which  also  determines  a  trial  value  of  M  from  which  in 
r  r 

turn  1^  is  computed.  Then  by  means. of  a  "corrector”  formula  an 
improved  value  of  1^  is  computed.  This  process  is  repeated  iteratively 
until  the  improved  values  is  identical  with  the  trial  value.  The  first 
trial  values  are  estimated  by  means  of  "predictor**  formulas. 


"Predictors" 


I2 

+  v 

(224) 

X3  = 

*2 

+  v 

§x 

(225) 
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(227) 
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X(i  +  1)  =  I(i-1)  +  3 


V f.  ,  +  41*.  +  I*,,  .v 
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(229) 
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In  computing  the  value  of  M  for  a  given  value  of  I,  the 
correction  to  be  added  to  a  trial  value  of  M  is 


§M 


M(1  +  v*^)(l  +  *jr  M2) 
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where  T  and  M  are  known  functions  of  x. 


The  results  of  Articles  11  and  12  to  be  tabulated  are, 
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M,  P 

i  l 


(234) 
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Plots  of  the  following  functions 


1  +  vM 


M  ^  1  M2 


I**!-,  M  (1  +  YM2)(1  +  *jr  M2), 


1  +  yM 


(235) 


[l  +  K2] 


useful  in  performing  the  computations  according  to  Articles 


12.11  and  12.12. 


12.13 


At  x  =  x3  identify 


y3J  ^  “  %  ,  *T?h  *  Ps3  ’  M3 
from  the  results  of  Articles  12.10,  12.11,  12.12  compute: 
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(239) 


F  =  2qv 
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In  241  ,  the  value  of  occurs  at  x  *  x^  and  the  values 


Ff,  2  occur  at  the  two  previous  stations  x. 


12.14  Determine  1^,  1^,  1^  and  from  the  system  of 


equations , 
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where 
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r2 


-  (Y  +  1)  Mx4 
L4(M12  -  1) 


1 


i 


3 


(247) 

(248) 

(249) 

(250) 

(251) 


In  order  to  solve  the  system  of  equations  242,  243,  244,  245,  and 
246,  compute  (y^  -  y^/l^  and  substitute  in  which  determines 

\x  .  Then,  compute  (y^  -  y^)/^  from  244  and  hence  (y^  - 

is  known  from  242  .  From  these  results  compute  l^j  l^j 

In  computing  the  value  of  {y^  -  y^)/l^  from  244  the  correction 
to  be  added  to  a  trial  value  of  (y^  -  y^/l^  is 
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(252) 


Finally,  compute 


3*c  =  (12  +  I3  +■  14)  _  _ ^ _ o  C  , 
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The 

12.13  are  to 


The 


'S’  V'  % 


computations  of  Articles  12.10,  12.11,  12.12,  and 

be  carried  to  the  x  *  x~  at  which  5x  =0. 

3  c 

input  data  for  Articles  12.13  and  12,14  include  besides 
*,  also  the  constants 


q, 


1 


The  results  of  Article  12.14  to  be  tabulated  are: 


y3>  y4>  l2>  V  V  V 

all  at  the  value  of  for  which 
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12.15  Compute,  finally 
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13.  SOLAR  ENERGY  COLLECTION1 

A  major  part  of  the  effort  on  this  contract  has  been  concerned  with 
investigating  lightweight,  high  efficiency  solar  energy  collectors,  which 
can  be  combined  with  a  suitable  heat  exchanger  for  the  purpose  of  absorbing 
energy  from  the  sun  and  transferring  it  to  a  working  fluid  with  the  highest 
possible  efficiency.  A  great  deal  of  research  effort  has  been  expended  for 
many  years  in  the  development  of  various  types  of  solar  energy  collectors. 
However,  no  devices  had  previously  been  developed  with  performance  approaching 
the  requirements  of  this  program.  The  solar  energy  collectors  to  be  used  for 
the  present  application  must  conform  to  extremely  rigid  requirements  if  the 
system  is  to  have  any  chance  of  successful  operation.  They  must  be  light¬ 
weight,  operate  with  minimum  radiation  losses,  cause  minimum  disturbance 
of  the  air  flow,  and  be  able  to  withstand  the  expected  environmental  con¬ 
ditions.  The  types  of  collector  now  being  considered  appear  to  have  the 
potential  for  accomplishing  these  objectives. 

13 . 1  General  Characteristics  of  Solar  Energy  Collectors 

An  irradiated  surface  tends  to  become  heated  when  its  optical 
absorption  is  high,  but,  upon  being  heated,  will  tend  to  cool  because  of 
thermal  reradiation.  The  spectral  distribution  of  the  absorbed  and  the 
re-emitted  radiation  depends  on  the  temperature  and  optical  absorptivity  of 
the  source  of  each.  Thus,  according  to  Planck* s  Law,  a  surface  at  temper¬ 
ature,  T,  emits  energy  with  intensity 

IW  -  «(V  (« czM-i yl  ^  c 

cm 

of  radiation  at  wavelength  \,  where  and  are  constants.  The  proportion 
of  incident  energy  at  wavelength  K  which  a  surface  absorbs  and  converts  into 


1  The  analyses  presented  in  Section  4  form  a  self-contained  unit  and  do  not 
necessarily  conform  with  the  nomenclature  given  in  the  front  of  this  report. 
The  symbols  used  in  this  section  arc  explained  as  they  are  introduced. 
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heat  is  given  by  its  absorptivity ,  a(/\.).  Thfe  only  other  physical  principle 
needed  to  formulate  the  problem  is  Kirchhoffts  Law,  which  states  that  for 
any  surface  whatsoever,  the  absorptivity  equals  the  emissivity  at  a  given 
wavelength. 

a(A.)  =  e(A.)  (257) 

Therefore,  the  first  problem  is  to  state  the  optimum  absorptivity  function, 

a  (a.  )  . 

Any  practical  solar  absorber  not  involving  direct  concentration 
of  the  sunrs  rays  will  attain  a  maximum  temperature  only  a  fraction  that  of 
the  sun.  The  spectral  energy  distribution  of  the  absorbed  radiation  will  be 
largely  distinct  from  that  of  the  reemitted  radiation  since  the  solar  spectrum 
is  approximately  that  of  a  black  body  at  550G°C,  with  a  peak  at  about  a.  «  0,5 
microns,  while  the  peak  reradiation  from  the  absorber  occurs  at  a  wavelength 
of  a  few  microns.  If  these  two  spectra  did  not  overlap  at  all,  then  the 
optimum  a  (A.)  would  obviously  be 

a  (A.)  *  1  in  the  solar  region,  and 
a  (A.)  -  0  in  the  reradiation  region; 

so  that  all  of  the  incident  energy  would  be  absorbed  and  none  reemitted.  As 
it  happens,  there  is  some  overlap  as  shown  in  Fig.  28a.  In  order  to  retain 
the  maximum  energy,  an  absorber  with  a  "cutoff"  characteristic  like  that 
shown  in  Fig.  28a  is  desirable. 

The  problem  is  now  to  find  the  best  value  for  the  cutoff  wavelength, 

A.  ,  which  depends  in  turn  on  the  resulting  absorber  temperature,  T.  If  the 

c 

"trapping  efficiency",  f,  is  defined  as  the  ratio  between  the  incident  energy 
converted  into  heat  and  the  total  incident  energy  (in  the  equilibrium  situ¬ 
ation)  then, 


where 


k^ 

[  isc^>^ 


^  is(\)dx 


(259) 


is  the  fraction  of  solar  radiation  absorbed,  and 

oo 
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iT0,)d(k) 


C  IT(>.)dV 


(260) 


is  the  fraction  of  black  body  reradiation  which  is  not  allowed  to  escape, 
IC(X)  is  the  energy  density  in  the  solar  spectrum  and  I  (A.)  is  that  in  the 
reradiation  spectrum.  A  plot  of  trapping  efficiency  vs.  temperature  is  given 
in  Fig.  29.  The  optimum  value  of  k^  is  then  found  by  selecting  the  desired 
temperature  curve  and  noting  its  maximal  point.  For  the  temperature  range 

shown  here,  k  should  thus  be  about  2.5  microns. 

9  c 

No  actual  absorber  can  attain  the  characteristics  shown  in 
Fig.  28a.  A  more  realistic ,  although  still  somewhat  idealized ,  collector 
characteristic  curve  is  shown  in  Fig.  28b.  The  trapping  efficiency,  f,  is 
defined  as  before: 


However,  f^  and  f^ 


are  now: 


als(\)d\ 


(K)dk 


1  oo 
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A. 

rc 


\  (1  -  €>IT(\)dA.  +  \  (1  -  €*)lTt/,)dA. 


CD 
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(262) 


The  quantities  a,  a* ,  €  and  €T  are  defined  in  Figure  28b.  The  optimum 
value  of  X^  is  found  as  before  for  any  particular  values  of  a(X)  and  e(X). 
Although  no  actual  collector  would  have  the  precise  "high-pass"  character- 

J 

is tic  illustrated  in  Fig.  28b,  this  serves  to  demonstrate  the  type  of  filter 
behavior  desired  for  solar  energy  collection  purposes.  An  kbsorber  becomes 
more  nearly  ideal  as 


a,  €->1.0  for  X  <  X  and  a?  ,  €r  — *  0  for  X  >  X 

For  a  real  surface  with  oi  and  €  continuously  varying  functions  of  X,  it  is 
necessary  to  integrate  the  products  a(X)Ig(X)  and  Q.  -  e(X)j  ^(X)  over  the 
entire  wavelength  range. 

There  are  many  types  of  surfaces  which  exhibit  the  absorption- 
emission  characteristics  illustrated  in  Fig.  28.  Among  these  are: 

1.  Metallic  smoke  deposits 

2.  Chemically  produced  surface  coatings  on  reflecting  metal 
subsurfaces 

3.  Interference  filters 

Many  examples  of  Types  1  and  2  have  been  produced  in  the  Electro -Optical 
Systems  laboratories  during  the  present  program.  Methods  of  producing  these 
surfaces,  as  well, as  performance  results,  are  discussed  in  the  succeeding 
sections . 

Although  some  have  very  promising  characteristics,  the  severe 
requirements  of  the  present  application  may  make  it  necessary  to  augment 
their  performance  by  making  use  of  the  fact  that  the  incident  solar  radiation 
is  very  nearly  parallel,  while  the  reradiation  from  the  absorbing  surfaces  is 
diffuse.  Since  the  incident  radiation  is  essentially  parallel,  it  can  be 
focussed  through  a  small  aperture  before  being  absorbed  by  a  collecting 
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surface.  The  diffuse  reradiation  then  "sees"  only  a  small  area  from  which 
to  escape.  Devices  which  operate  in  the  above  manner  are  termed  geometric 
radiation  traps.  Examples  of  these,  have  been  built  and  tested  at  Electro- 
Optical  Systems.  In  order  to  gain  the  maximum  performance,  the  geometric 
trap  is  used  in  conjunction  with  a  spectrally  selective  absorbing  surface. 
Many  different  configurations  are  possible  using  the  basic  focussing 
principle  to  reduce  reradiation. 

Absorbing  surfaces  using  the  intereference  principle  may  have 
very  good  absorbing-emitting  characteristics  when  used  alone  with  normal 
incidence  light.  However,  it  is  believed  by  the  ElectrorOptical  Systems* 
staff  that  their  usefulness  in  conjunction  with  a  geometric  radiation  trap 
would  be  limited.  Therefore,  the  emphasis  during  the  present  program  has 
been  based  on  other  types  of  spectrally  selective  absorbing  surfaces.  The 
general  operating  characteristics  of  the  various  types  of  absorbing  surfaces 
and  geometric  traps  are  discussed  in  the  following  sections. 

13 . 2  Spectrally  Selective  Surfaces 

Solar  energy  research  in  recent  years  has  resulted  in  several 
types  of  surface  coatings  with  the  absorptivity  versus  wavelength  character¬ 
istics  discussed  in  Section  13.1.  These  characteristics  include  high 
absorptivity  in  the  visible,  low  emissivity  in  the  infrared,  and  a  suitable 
cutoff  wavelength  which  maximizes  total  energy  absorption. 

Aside  from  optical  interference  filters  (discussed  in  Section 
13.4)  these  surfaces  are  prepared  by  metallic  or  non-metallic  smoking 
processes  and  chemical  or  electrochemical  deposition  processes. 

When  embarking  on  an  experimental  program  in  spectral  selectivity, 
it  was  Electro-Optical  Systems*  aim  to  duplicate  and  verify  the  results  of 
other  workers  in  the  field  as  V7ell  as  to  develop  new  coatings  which  would 
show  improved  characteristics.  Both  these  aims  have  been  accomplished  on  the 
present  program.  This  section  includes  a  detailed  discussion  of  the  methods 
of  preparation,  general  physical  characteristics,  and  performance  of  the 
coatings  prepared  in  the  EOS  laboratory. 
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Several  types  of  surfaces  with  suitable  cutoff  wavelengths  have 
been  investigated.  One  or  more  processes  in  each  of  the  following  areas  have 
been  tried  and  tested: 

1*  Gold  smoke  surfaces 

2.  Oxidized  sprayed  coatings. 

3.  Electrochemically  deposited  coatings 

4.  Metal  oxide  coatings  obtained  by  oxidizing  the  substrate 
metal 

5.  Non-metallic  smoked  coatings 

6.  Metal  oxide  coatings  on  glass  (NESA  glass) 

7.  Chemically  plated  coatings 

These  types  of  coatings  or  surfaces  are  described  in  greater  detail 
below.  The  surfaces  are  evaluated  on  the  basis  of  reflectance  and  transmit¬ 
tance  measurements  using  a  13U  Perkin-Elmer  spectrophotometer.  All  samples 
were  evaluated  from  below  0.3  to  above  14  microns. 

Figure  31  summarizes  the  results  obtained  in  the  visible  and 
infrared  for  the  best  black  surfaces.  Low  reflectance  in  the  visible  (high 
absorptance  or  transmittance)  and  high  reflectance  in  the  infrared  (low 
emittance  or  transmittance)  are  the  spectral  characteristics  desired.  An 
inflection  point  between  1.5  -  3.5  microns  is  optimum  (see  Fig.  31). 

Figures  32  and  33  summarize  the  reflectance  and  transmittance  versus  wave¬ 
length  data  for  the  best  metal-oxide  coated  glass. 

The  reflectance  measurements  are  specular  reflection  only  (total 
reflection  equals  specular  plus  diffuse)  from  an  incidence  angle  10°  off 
the  normal.  It  is  not  anticipated  that  the  wavelength  dependence  of  reflec¬ 
tion  is  a  strong  function  of  the  type  of  reflection  (specular  normal  instead 
of  a  total  hemispherical  reflection).  However,  since  percent  specular 
reflection  is  a  function  of  surface  roughness  and  angle  of  incidence,  it 
would  be  desirable  to  obtain  total  hemispherical  reflection  data.  However,* 
the  equilibrium  temperature  tests  described  in  Section  13.3  serve  as  an 
over-all  evaluation  of  total  hemispherical  reflection.  The  reflection 
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data  obtained  from  the  spectrophotometer,  therefore,  serve  as  adequate 
first  approximations  to  the  nature  of  the  surface. 

It  should  be  emphasized  that  the  spectral  properties  of  many 
of  the  surfaces  are  strongly  dependent  on  coating  thickness.  The  commer¬ 
cial  processes  investigated  were  generally  prepared  in  accordance  with 
manufacturer^  directions.  The  experimental  controls  were  then  varied  to 
produce  surfaces  with  optimum  spectral  characteristics.  During  this  work 
many  surfaces  were  improved  by  thinning  the  surface  film. 

13.2.1  Theoretical  Discussion 

Several  years  ago  it  was  pointed  out  by  Tabor  (Ref.  12) 
that  the  efficiency  of  the  collecting  surfaces  used  as  solar  heating  devices 
could  be  increased  if  they  were  coated  with  a  material  which  absorbs-  in  the 
range  of  the  solar  spectrum,  but  is  transparent  at  long  wavelengths,  corres¬ 
ponding  to  the  temperature  to  which  the  collector  is  heated  by  the  sun.  If 
this  coating  were  backed  by  a  metallic  reflector  of  low  absorptivity  (and 
emissivity)  in  the  infrared,  then  incoming  solar  radiation  would  be 
absorbed  by  the  coating  and  the  absorbed  energy  transferred  to  the  metallic 
substrate  by  thermal  conduction  in  the  solid  state.  Due  to  the  poor 
emissivity  of  the  metal  surface,  very  little  of  this  energy  is  reradiated 
from  the  metal  at  the  temperature  (usually  a  few  hundred  degrees)  to  which 
the  surface  is  raised.  Thus  radiation  losses  from  the  collector  surfaces 
would  be  almost  eliminated,  metal  emissivities  being  only  3-4  percent. 

The  surfaces  used  by  Tabor  and  subsequent  workers  in 
this  field  were  black  metallic  oxides,  usually  deposited  on  aluminum 
backing.  Copper  oxide  (CuO)  has  been  investigated  rather  thoroughly  by 
Unger  in  his  M.I.T.  thesis  (Ref.  13)  and  further  investigations  have  been 
carried  out  at  EOS  under  the  present  contract. 

The  work  done  so  far  has  been  limited  to  the  development 
of  satisfactory  chemical  methods  for  depositing  oxides  or  other  compounds 
on  the  metallic  surfaces,  and  the  measurements  of  reflectivities  and 
absorptivities  of  the  surfaces  in  the  visible  and  near  infrared  ranges. 
However,  Unger  has  mentioned  the  fact  that  optical  absorption  in 


74 


semiconductor  materials  is  due  to  an  excitation  of  electrons  from  the  valence 
band  into  the  conduction  band,  across  the  energy  gap  existing  in  the  semi¬ 
conductor.  The  absorption  edge  is  given  by  radiation  whose  photon  energy 
is  equal  to  the  minimum  energy  gap  between  valence  and  conduction  bands. 

If  surface  scattering  effects  are  ignored,  then  radiation  of  energy  above 
the  "gap  value"  will  interact  with  electrons  excited  into  the  conduction 
band  and  absorption  will  occur  as  a  result  of  transitions  in  this  band. 

Photons  of  energy  less  than  the  'energy  gap'  will  not 
excite  electrons  into  the  conduction  band  and  the  material  will  be  trans¬ 
parent  to  them.  Since  E  =  h,  the  low  energy  photons  correspond  to  low 
frequency  (long  wavelength)  radiation.  This  is  the  qualitative  explanation, 
in  terms  of  solid  state  properties,  of  the  spectrally  selective  behavior  of 
certain  surfaces.  This  simple  method  is  complicated  by  two  factors--f irst , 
the  effects  of  impurities  in  the  semiconducting  surface;  secondly,  the 
scattering  action  of  geometrical  irregularities  in  the  surface. 

In  terms  of  the  band  theory  of  semiconductors,  the  first 
feature  was  treated  by  the  transistor  developers,  Bardeen,  Shockley  and 
associates  (Ref.  14).  Atoms  of  a  material  with  a  valence  one  greater  or  one 
less  than  that  of  the  pure  semiconductor  are  added  to  it,  producing  a 
surplus  or  deficit  of  electrons,  respectively,  in  the  lattice  structure. 

In  terms  of  the  band  theory,  the  impurities  introduce  electrons  whose  energy 
levels  lie  in  the  gap  between  valence  and  conduction  bands.  Under  the 
influence  of  absorbed  photons  of  appropriate  energy,  electrons  may  be 
excited  from  these  'impurity1  levels  into  the  conduction  band  or  from  the 
valence  band  to  an  impurity  level,  leaving  a  'hole'  in  the  valence  band 
which  migrates  as  a  current  carrier  in  the  same  way  as  an  electron  in  the 
conduction  band. 

Clearly,  the  presence  of  these  impurity  levels  in  the 
energy  gap  can  greatly  modify  the  absorption  characteristics  of  a  semiconductor, 
since  the  effective  width  of  the  energy  gap  is  changed  (generally  reduced)  . 

Such  an  effect  appears  to  occur  in  the  ARL  supported  experiments  on  CdS  as 
a  photovoltaic  cell  material.  However,  no  systematic  investigation  of 
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the  effects  of  such  impurities  on  absorption  properties  of  semiconductors 
and  the  relation  of  these  properties  to  their  usefulness  as  solar  collectors, 
has  yet  been  carried  out. 

The  second  complicating  factor,  mentioned  above,  is  the 
light  scattering  effect  due  to  small  surface  irregularities.  It  is  con¬ 
venient  to  summarize  the  relevant  work  on  this  topic.  The  problem  of  the 
scattering  of  electromagnetic  waves  by  a  sphere  of  arbitrary  size  and  „ 
dielectric  constant  was  solved  in  1908  by  Mie  (Ref.  15),  using  the  standard 
methods  of  electromagnetic  field  theory.  His  results  give  intensity  of 
scattered  light  in  terms  of  initial  intensity,  as  a  cumbersome  and  slowly 
convergent  series  expression.  However,  numerical  values  of  this  have  been 
computed  over  a  wide  range  of  parameters  and  some  of  the  results  are 
summarized  in  van  de  Hulst’s  book  (ref.  16). 

It  rs  round  that ,  as  a  result  of  interference  during  the 

scattering  process,  a  certain  percentage  of  light  is  extinguished.  Final 

numerical  results  appear  as  graphs  of  the  extinction  percentage,  Q,  against 

2jt  r 

a  parameter  x  of  value  — r - ,  where  r  is  the  radius  of  the  scattering 

X 

sphere,  and  X  is  the  wavelength  of  the  incident  light. 

Graphs  of  Q  against  x  have  been  made  for  both  dielectric 
and  conducting  spheres.  They  indicate  rather  sharp  peaks  in  the  extinction 
coefficient  at  a  certain  value  of  x.  This  critical  value  differs  considerably 
for  conducting  and  non-conducting  spheres,  being  considerably  smaller  in  the 
case  of  the  former.  In  all  cases,  the  extinction  coefficient  becomes  zero 
rapidly  as  x  goes  to  zero.  Beyond  the  critical  values  of  x,  the  extinction 
curve  drags  somewhat  and  flattens  out.  The  results  of  these  scattering 
calculations  may  be  applied  to  a  rough  surface  by  setting  up  an  idealized 
model  of  this  surface,  consisting  of  a  perfectly  flat  backing  dotted  with 
spheres  forming  a  partially  opaque  optical  layer  immediately  in  front  of  the 
backing.  On  an  actual  surface,  the  various  irregularities  will  not  be 
spherical,  nor  will  they  necessarily  be  independent  of  their  neighbors, 
so  far  as  scattering  effects  go.  Hence  the  discussion  here  can  only  be  con¬ 
sidered  a  zero- order  approximation. 
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This  result  follows.  When  the  ratio  r/X  is  in  the  range 
0.1  to  1.0,  the  extinction  coefficient  is  at  a  maximum.  When  r/X  is  less 
than  this,  a  rather  sharp  cutoff  in  extinction  should  be  noted  (in  the 
longer  wavelengths  for  constant  r).  When  r/X  becomes  greater  than  this, 
one  can  expect  a  moderate  increase  in  transmission  of  light.  It  may  be 
noted  that  the  optical  properties  of  the  colloidal  gold  deposit  studied  by 
Harris  (Ref.  17)  (and  currently  at  EOS)  are  in  agreement  with  this  statement. 
For  the  limit  where  X<$^  r  one  gets  scattering  from  particles  in  accordance 
with  ordinary  rules  of  geometrical  optics.  By  multiple  reflection  from  its 
surface  irregularities,  the  normal  absorptivity  of  the  medium  is  increased. 
This  seems  to  explain  the  effect  of  surface  roughening  in  increasing 
absorptivities  (and  hence  emissivities)  of  surfaces. 

It  is  probable  that  the  radiation  reflectivity  of  certain 
white  oxides  and  ceramic  materials  is  due  to  scattering  effects.  As  x 
decreases  beyond  the  critical  value  (or  as  X  increases),  one  *'ould  get  an 
increase  in  the  extinction  coefficient  near  a  resonant  peak  and  for  certain 
values  of  the  refractive  index,  this  peak  can  be  broad  enough  to  cover  a 
considerable  wavelength  region.  As  far  as  is  known,  no  systematic  investi- 
fation  of  the  relation  between  surface  reflectivities  as  various  wavelengths 
and  the  size  of  the  diffracting  obstacles,  namely  the  irregularities  in  the 
surface,  has  been  undertaken.  Such  an  investigation,  together  with  a  quanti 
tative  theoretical  discussion,  would  be  needed  to  settle  the  problem. 

Finally,  it  may  be  pointed  out  that  the  work  of  Tabor  and 
others  on  spectrally  selective  semiconductor  surfaces  has,  to  date,  been 
carried  out  at  low  surface  emission  temperatures,  which  are  appropriate  for 
extended  solar  collector  surfaces  under  atmospheric  conditions.  However, 
for  solar  power  sources  under  space  conditions,  a  considerably  higher 
concentration  of  solar  energy  is  necessary,  and  thus  temperatures  in  the 
range  2000°  to  3000°F  will  be  encountered.  The  problem  of  developing 
suitable  spectrally  selective  coatings  for  use  under  these  conditions  does 
not  seem  to  have  received  much  attention  as  yet.  The  only  high  temperature 
semiconductor  whose  optical  properties  have  been  investigated  is  silicon 
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carbide.  Data  on  high  temperature  oxides  and  pyrites -type  compounds  are  not 
presently  available. 

Under  high  temperature  conditions,,  moreover,  a  significant 
portion  of  absorbed  energy  may  go  into  activating  the  thermionic  emission 
of  electrons.  How  far  the  thermionic  emission  characteristics  may  be 
affected  by  the  photoactivated  transfer  of  electrons  from  valence  to  con¬ 
duction  bands,  and  the  effect  of  impurities  in  this  process,  is  not  known. 
The  mechanism  of  conversion  of  radiation  energy  to  energy  of  thermionically 
emitted  electrons  requires  investigation  as  this  process  is  part  of  the 
mechanism  involved  in  optical  absorption  by  semiconductors. 

Summarizing,  this  discussion  has  outlined  two  distinct 
spectral  selectivity  effects.  One  is  due  to  scattering  by  surface 
irregularities,  the  other,  and  more  interesting,  to  the  peculiar  electronic 
structure  of  semiconductors.  The  most  promising  line  of  further  work  seems 
to  lie  in  a  study  of  this  latter  effect  at  higher  temperatures,  together 
with  the  closely  related  thermionic  emission  properties  of  an  absorbing 
semiconductor . 

13.2.2  Description  and  Preparation  of  Surfaces 

Various  references  in  the  literature  (Refs.  13  and  18 
through  26)  have  been  investigated  with  the  intention  of  duplicating  the 
state-of-the-art  and,  if  possible,  of  making  improvements  in  spectrally 
selective  surfaces.  These  include  black  or  near  black  coatings  on  aluminum, 
copper,  magnesium,  steel,  platinum,  nickel,  and  glass;  as  well  as  electri¬ 
cally  conductive  coatings  on  glass. 

This  section  describes  the  physical  properties  and  methods 
of  preparation  of  these  surfaces.  Following  sections  discuss  their  spectral 
properties  and  suitability  for  use  in  solar  energy  collectors. 

13.2.2.1  Gold  Smoke  on  Glass 

The  gold  smoke  filter  derives  its  spectral 
properties  from  the  mechanical  structure  of  its  surface  which  is  composed 
of  many  microscopic  gold  droplets.  Differential  scattering  of  short  and 
long  wave  radiation  gives  the  characteristic  of  high  absorptance  of  solar 
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radiation  but  low  emittance  of  infrared.  The  properties  and  spectral 
characteristics  of  gold  smokes  are  described  in  Refs.  17,  18,  and  19. 

Experimental  work  has  been  undertaken  at 
Electro-Optical  Systems  which  is  aimed  at  producing  improved  forms  of  the 
gold  smoke  filter.  The  gold  smoke  surfaces  produced  during  this  program 
have  been  deposited  on  a  glass  substrate  in  order  to  facilitate  fabrication. 
However,  if  they  were  to  be  used  in  a  collector  they  would  be  put  on  an 
aluminum  or  aluminized  glass  base  in  order  to  take  advantage  of  the  low 
metallic  emissivity  in  the  infrared.  Gold  smoke  filters  are  prepared  in 
the  following  manner. 

Gold  wire  was  fused  to  a  molybdenum  boat  in 
vacuum  deposition  system  at  a  pressure  of  5  microns  and  at  white  heat  in 
each  case.  A  2"  x  2"  glass  sheet,  which  in  some  cases  was  preheated  in  an 
air  oven,  was  then  placed  at  3M  or  3  l/2n  above  the  boat,  and  the  system 
evacuated  until  an  air  pressure  of  between  80  and  95  microns  was  measured. 
Pure  dry  nitrogen  gas  was  then  admitted  through  a  needle  valve  until  a 
total  pressure  of  2  mm  of  Hg  was  measured.  The  boat  was  then  electrically 
heated  by  specific  currents  for  varying  amounts  of  time.  The  system  was  then 
shut  down,  the  glass  sheet  removed  and  its  transmittance  measured  in  the 
range  from  0.4p  to  about  4.0fi.  Among  the  25  samples  thus  prepared,  five 
promising  gold  smoke  deposits  were  found.  A  tabulation  of  the  variables  for 
these  five  samples  is  shown  below,  while  their  transmission  spectra  are 
reproduced  in  Fig.  34. 


Sample 

Number 

Air  Pressure 
(microns) 

Boat  to 
Glass 
Distance 

Filament  Setting 
(arbitrary 
units) 

Deposition  Time 

Glass  Temp. 
(Preheating) 

5 

90 

3” 

15 

40  sec 

20°C 

6 

85 

.3" 

15 

40  sec 

20°C 

19 

85 

3  1/2" 

15 

120  sec 

20°C 

20 

85 

n  1  /  >-»  II 
i./  i. 

1  c 

I.J 

until  light  gray 

t 

i.  L.\J  O 

23 

85 

3  1/2" 

15 

until  dark  gray 

120°C 
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One  of  the  foremost  problems  in  investigating 
gold  smoke  deposits  is  that  of  reproducibility.  Several  samples  were 
prepared  under  apparently  identical  conditions,  but  showed  widely  differing 
spectral  characteristics.  More  conclusive  results  will  probably  be  obtained 
by  the  use  of  boats  of  smaller  size  and  greater  depth  so  that  during  the 
fusing  stage  of  the  process  the  gold  will  adhere  to  the  boat  in  a  more 
uniform  and  reproducible  pattern. 

One  aspect  of  this  work  has  been  particularly 
promising,  namely  the  temperature  stability  of  the  deposits  *  They  have 
maintained  their  spectral  characteristics  and  physical  integrity  at  a 
temperature  of  140°C,  whereas  previous  deposits  prepared  by  other  workers 
were  found  to  disintegrate  at  120°C. 

13 . 2 . 2 . 2  Sprayed  Cupric  Oxide  on  Aluminum 

These  surfaces,  first  prepared  by  Unger  (Ref. 

13),  are  cupric  oxide  coatings  formed  by  spraying  cupric  nitrate  solution 
onto  heated  samples  of  reflective  aluminum. 

Figure  31  indicates  that  Unger  surfaces  are 
among  the  best  spectrally  selective  surfaces  which  Electro-Optic^x  Systems 
has  found. 

The  optimization  of  high  solar  absorptivity  and 
low  infrared  emissivity  can  be  closely  controlled. 

Figure  35  from  Unger* s  thesis  indicates  the 
range  in  the  nature  of  surfaces  which  can  be  made,  expressed  in  integrated 
solar  absorptivities  (a)  and  integrated  80°C  black-body  emissivities ,  (e). 
Additional  analysis  is  required  to  determine  the  optimum  obtainable  a  +  g 
for  any  particular  requirement,  since  e  increases  as  a  increases. 

The  Unger  surface  is  a  cupric  oxide  on  aluminum 
coating.  It  is  strongly  adherent,  stable  at  high  temperatures,  and  can  be 
prepared  under  very  reproducible  conditions.  A  cupric  nitrate  solution  is 
sprayed  onto  a  heated  aluminum  surface.  The  surface  is  then  baked,  converting 
the  cupric  nitrate  to  cupric  oxide. 
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The  spray  head  is  a  125  M  round  spray  head 
from  the  Spray  Engineering  Company,  100  Cambridge  Street,  Burlington, 
Massachusetts.  The  fluid  nozzle  is  designated  F2  and  the  air  nozzle 
F2081.  The  nozzles  are  integral  parts  of  the  spray  head.  The  aluminum  is 
99.99  percent  pure  (0  or  1/2  hard  temper)  obtained  from  the  Alcoa  Research 
Laboratories  in  New  Kensington,  Pa.  Unger  found  that  this  high  purity 
aluminum  was  necessary  in  order  to  insure  proper  surface  reflectivity  of 
the  aluminum  base.  The  aluminum  was  used  as  received,  requiring  no  surface 
treatment.  However,  the  four-inch  squares  which  were  used  in  the  heat  equi¬ 
librium  tests  had  to  be  drilled,  and  in  the  process  they  picked  up  oil  which 
was  cleaned  off  by  successive  dips  in  xylene,  redistilled  acetone,  and  ethyl 
alcohol . 

A  compressor  (with  a  filter  to  remove  oil  from 
the  airstream)  provides  the  airflow  at  5  -  20  psig.  The  aluminum  sample  is 
fastened  to  a  hot  plate  (a  thermocouple  is  also  fastened  to  the  hot  plate 
just  at  the  edge  of  the  aluminum  sample).  The  aluminum  is  placed  in  a 
vertical  position  for  horizontal  spraying  22  cm  from  the  nozzle  head.  The 
solution  is  fed  to  the  nozzle  from  a  constant  head  flask  placed  30  cm  above 
the  level  of  the  nozzle. 

Eight  variables  affect  the  low  temperature 
emission  of  the  Unger  surface  (absorptivity  for  sunlight  is  relatively 
insensitive  to  the  variables).  The  following  results  were  obtained  by 
Unger  and  were  generally  verified  by  this  laboratory. 

1 .  Concentration  of  Cupric  Nitrate  in  the  Spray 
0.0025  molal  cupric  nitrate  is  a  near  opti¬ 
mum  concentration  to  give  the  proper  particle 
size  of  deposit;  0.005  molal  is  acceptable. 

2.  Spray  Particle  Droplet  Size 

The  emissivity  increases  as  the  droplets 
decrease  in  size  from  60  to  12  microns.  This 
variable  is  controlled  by  the  spray  rate 
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discussed  below. 


3 .  Spraying  Rate 

The  optimum  spray  rate  is  about  150  cc/hour; 
the  acceptable  limits"  being  50  -  250  cc/hour. 
At  this  rate  the  droplet  size  and  spray  rate 
effect  on  the  drying  of  the  film  are  optimum. 
This  rate  is  obtained  by  the  screw  adjust¬ 
ment  of  the  fluid  nozzle. 

4.  Drying  Plate  Temperature 

The  hot  plate  can  be  kept  at  any  tempera¬ 
ture  high  enough  to  prevent  visible  wetting 
of  the  aluminum  surface.  Temperatures  should 
not  exceed  170°C,  since  the  nitrate  converts 
to  the  oxide  above  this  temperature .  The 
range  studied  was  from  80°  to  170°C. 

5.  Baking  Temperature 

Minimum  emissivity  and  maximum  solar  absorp¬ 
tivity  are  obtained  by  baking  at  350°C. 

At  230°C  on  one  extreme,  and  450°C  on  the 
other,  definite  detrimental  factors  set  in. 

6.  Time  of  Baking 

Samples  baked  for  at  least  twelve  hours 
had  lower  emissivities  than  samples  baked 
for  shorter  times.  No  effect  on  optical 
properties  of  the  surface  was  evidenced  by 
varying  the  rate  of  heating  or  cooling. 

7 .  Interval  Between  Completion  of  Spraying 
and  Baking 

The  adhesion  of  the  coating  to  the  aluminum 
is  good,  but  the  cohesion  of  the  cupric 
oxide  film  to  itself  is  adversely  affected 
by  a  lapse  in  time  before  baking.  If  baked 
immediately  after  spraying  the  cupric  oxide 
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film  will  resist  strong  vibrations  and 
air  currents.  If  alternately  baked  and 
sprayed,  the  coating  will  resist  vigorous 
massage. 

8.  Amount  of  Deposit  Per  Unit  Area 

2 

The  optimum  range  is  0.275  to  0.350  mg/ cm 

of  CuO.  For  an  optimum  spray  rate,  this 

indicates  a  spraying  time  of  about  one 

hour.  This  calculation  is  based  on  the 

fact  that  150  cc/hr.  spray  rate  gives 
2 

1 . 5  cc/ cm  hr . 

The  spray  covered  about  a  square  inch  of 
surface  uniformly.  In  order  to  cover  4M  x  4n  aluminum  squares  with  an 
Unger  coating  the  hot  plate  was  set  on  a  track  and  slowly  pulled  in  one 
direction.  Whereas,  the  spray  was  attached  to  a  cam  shaft  operated  off  a 
heart-shaped  cam  which  provided  a  constant  velocity  (except  at  the  ends  of 
the  stroke)  reciprocating  motion.  A  picture  of  this  experimental  apparatus 
is  given  in  Fig,  30. 

13 . 2 . 2 . 3  Cupric  Oxide  on  Copper  by  the  Ebonol  Process 
(Ref.  20) 

A  deep  black  surface  with  a  velvet-like 
appearance  can  be  obtained  on  copper  by  this  process.  The  procedure  for 
obtaining  this  surface  is: 

1.  Preparation  of  copper  base  (copper  used  is  99+ 

percent  pure): 

a.  Scrub  with  Dutch  cleanser  to  remove  grease. 

b.  Etch  in  10  N .  nitric  acid  to  uniformly  activate 
the  surface. 

c.  Pickle  in  10  percent  sulfuric  acid. 

2.  Preparation  of  Ebonol  bath: 

a.  90.0  grams  of  Ebonol  "C"  Special  are  dissolved 
in  500  ml.  water. 

b.  The  bath  is  heated  to  99  -  102°C. 
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3.  Preparation  of  surface: 

a.  The  copper  is  immersed  in  the  Ebonol  bath  for 
8  -  10  minutes, 

b.  The  black  surface  is  then  rinsed  in  warm  water. 

c.  The  surface  is  air-dried. 

This  patented  product,  Ebonol,  was  given  a 
cursory  analysis  and  determined  to  be  an  alkaline  hypochlorite. 

The  deep  black  of  this  surface  is  damaged  by 
touching  it,  as  the  fine  protrusions  of  cupric  oxide  are  collapsed  under  the 
slightest  force. 

13.2.2.4  Black  Chromate  on  Magnesium  (Ref.  21) 

A  dark  brown  coating  can  be  obtained  on 
magnesium  by  this  process.  The  preparation  procedure  is; 

1.  Preparation  of  magnesium  base: 

a.  Scrub  with  Dutch  cleanser  to  remove  grease. 

b.  Pickel  in  60  percent  hydrofluoric  acid  at 
21-32°C. 

2.  Preparation  of  chromate  bath: 

a.  A  bath  of  the  following  composition  is  prepared: 

Ammonium  sulfate,  (NH^^SO^  30.0  g/1 

Sodium  dichr ornate,  Ite^Cr  20^2^0  30.0  g/1 

Ammonium  hydroxide,  NH^OH  2.6  ml/1 

(Sp.  gr.  0.880) 

b.  The  pH  of  the  bath  is  5. 6-6.8. 

c.  The  bath  is  maintained  at  99  -  100°C. 

3.  Preparation  of  surface: 

a.  The  magnesium  is  immersed  in  the  chromate  bath 
bath  for  30  -  40  minutes. 

b.  The  coated  surface  is  rinsed  several  times  in 
warm  water. 

c.  The  surface  is  air-  or  over-dried  at  93°C  for 
5  to  15  minutes. 
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A  chocolate  brown  surface  was  the  darkest 
surface  which  could  be  obtained  on  pure  magnesium  by  this  process.  Any 
attempt  to  produce  a  darker  surface  only  resulted  in  a  powdery  non-adherent 
coating. 

13.2.2.5  Black  Nickel-Iron  Oxide  on  Steel,  by  the 
"Nickel  Pentrate"  Process  (Ref.  22) 

A  blue  gray-black  adherent  surface  can  be 
obtained  on  steel  by  this  process.  The  procedure  for  obtaining  this  surface 
is  as  follows: 

1.  Preparation  of  steel  base  (magnetic  steel): 

a.  Scrub  with  Dutch  cleanser  to  remove  grease. 

b.  Etch  in  10  -  20  percent  sulfuric  acid. 

2.  Preparation  of  Nickel  Pentrate  bath: 

a.  180  g.  of  Nickel  Pentrate  are  dissolved  in 
250  ml.  water. 

b.  The  bath  is  maintained  at  126°C.  (The  solution 
boils  at  126°C  although  the  recommended  temper¬ 
ature  at  which  to  maintain  the  bath  is  143. 3°C. 
This  discrepancy  evidently  arises  from  the  fact 
that  the  boiling  point  is  dependent  upon  the 
salt  concentration). 

3.  Preparation  of  surface: 

a  The  steel  is  immersed  in  the  Pentrate  bath 
for  15  minutes . 

b.  The  surface  is  rinsed. 

c.  The  surface  is  air-dried. 

In  the  chemical  process  of  oxidation,  the 
nickel  becomes  an  integral  part  of  the  black  oxide  surface  imparted  to  the 
steel.  Spectrographic  analysis  and  microscopic  examination  show  the  presence 
of  the  nickel  contained  in  the  black  oxide  finish.  The  appearance  of  this 
surface  is  similar  to  that  of  the  blueing  given  to  steel  guns. 

13.2.2.6  Black  Nickel  on  Copper  or  Aluminum  (Ref.  23) 
Blue-black  and  brass-colored  surfaces  can  be 
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formed  on  copper  by  this  electrolytic  deposition  method.  The  preparation 
procedure  is  as  follows: 

1.  Preparation  of  copper  base: 
a.  Scrub  with  Dutch  cleanser  to  remove  grease. 

2.  Preparation  of  black  nickel  bath: 

Nickal  sulfate,  NiSO^ 

Sodium  Thiocyanate,  NaSCN 
Zinc  sulfate,  ZnSO^ 

Lead  acetate,  Pb(CoH^02)2 

3.  Preparation  of  surface: 

Operating  conditions: 

Temperature: 

Voltage: 

Current  density: 

Anodes: 

Agitation: 

A  gray-black  surface  can  be  obtained  on  aluminum 
by  this  method.  The  procedure  for  obtaining  the  surface  is  the  same  as  that 
for  black  nickel  on  copper.  In  order  to  increase  the  current  density,  which 
is  very  low  because  of  the  thin  oxide  film  which  makes  the  aluminum  cathode 
passive,  the  aluminum  is  given  a  zincate  dip  which  deposits  a  thin,  uniform 
coat  of  zinc  on  the  aluminum.  The  current  density,  by  this  means,  was 
raised  to  one-half  the  recommended  current  density. 

13.2.2.7  Black  Nickel  on  Nickel  and  Nickel-Plated 

Copper 

Actual  tests  were  conducted  which  proved 
that  black  nickel  deposits  on  copper  surfaces  do  not  suffer  any  loss  in 
optical  properties  when  heated  to  170°C  in  a  300  micron  vacuum  test  chamber. 
Hov7ever,  the  black  nickel  on  nickel  system  has  been  studied  in  case  a  system 
with  exceptionally  high  resistance  to  thermal  deterioration  would  he  required. 
The  procedure  for  obtaining  the  black  nickel  on  nickel  surface  is  as  follows: 


97.5  g/1 

70.5  g/1 
45.0  g/1 

11.25  g/1 


27  C-38  C 
0.75-1.5  volts 
1-2  amps/ ft2 
nickel  or  carbon 
cathode,  when  necessary 
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1.  Preparation  of  nickel  base: 

a.  99+  percent  pure  nickel  is  polished  to  a  near 
mirror-bright  finish. 

b.  Scrub  with  Dutch  cleanser  to  remove  grease. 

(See  13.2.2.6) 

2.  Preparation  of  surface; 

(See  13*2.2.6) 

Black  nickel  was  deposited  both  on  nickel 
and  nickel-plated  copper  in  order  to  determine  if  there  were  any  differences 
in  emlssivity  characteristics  between  the  plate  nickel  and  the  electro- 
lytically  deposited  nickel.  Also  it  would  be  useful  to  determine  whether 
black  nickel  on  nickel  would  be  more  resistant  to  oxidation  than  black 
nickel  on  nickel-plated  copper,  since  it  is  well  known  that  electrolytically 
deposited  nickel  is  quite  porous  unless  special  techniques  of  intermittent 
plating  and  buffing  are  employed.  The  procedure  for  obtaining  black  nickel 
on  a  nickel-plated  copper  surface  is  as  follows: 

1.  Preparation  of  nickel-plated  copper: 

a.  99+  percent  pure  copper  is  polished  to  a 
mirror-bright  surface. 

b.  The  copper  is  given  a  flash  coat  of  nickel  in 
a  Watts* -type  nickel  plating  bath. 

2.  Preparation  of  black-nickel  bath: 

(See  13.2.2.6) 

3.  Preparation  of  surface: 

(See  13.2.2.6) 

13.2.2.8  Cupric  Oxide  on  Copper  by  the  "Dura-Black" 

Process  (Ref.  24) 

A  coating  similar  in  appearance  to  that 
produced  by  the  Ebonol  procedure  can  be  obtained  by  this  process,  which  is 
one  of  single  imnersion  in  an  alkaline,  organic  solution.  This  surface  does 
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not  have  the  same  depth  of  black.,  however,  the  minute  protrusions  of 
copper  oxide  are  the  same,  as  is  the  ea.se  with  which  they  can  be  damaged. 
The  basic  black  coating,  nevertheless,  is  extremely  hard  and  adherent  and 
will  withstand  bending,  twisting,  and  forming.  The  procedure  for  obtaining 
this  surface  is  as  follows; 

1.  Preparation  of  copper  base; 

a.  Scrub  with  Dutch  cleanser  to  remove  grease. 

b.  Bright  dip. 

2.  Preparation  of  Dura-Black  bath: 

a.  90.0  g.  of  Dura-Black  are  dissolved  in  1000 
ml.  water. 

b.  The  black  surface  is  then  rinsed  in  warm  water. 

3.  Preparation  of  surface: 

a.  The  copper  is  immersed  in  the  Dura-Black  bath 
for  4-5  minutes. 

b.  The  black  surface  is  then  rinsed  in  warm  water. 


13.2.2.9  Platinum  Black  on  Platinum  and  Nickel  (Ref.  25) 


A  velvet- like  black  coating  can  be  obtained 
on  platinum  by  standard  methods  of  electrolytic  deposition.  The  procedure 
for  obtaining  this  surface  is  as  follows: 

1.  Preparation  of  platinum  base: 

a.  Scrub  with  Dutch  cleanser  to  remove  grease. 

b.  Dip  in  aqua  regia. 

c.  Rinse  well. 

2.  Preparation  of  platinizing  solution: 

a.  The  solution  is  composed  of  3.0  percent  chloro- 
platinic  acid  and  0.2  percent  lead  acetate. 

3.  Preparation  of  surface: 


Operating  conditions: 
Temperature 
Voltage 
Anodes 
Agitation 


Room  temperature 
4.0  volts 

Carbon  or  platinum 
Cathode,  when  necessary 
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The  procedure  which  is  generally  recommended 
is  to  reverse  the  voltage  every  minute  until  a  thick  coating  is  obtained; 
however,  if  a  power  supply  with  an  AC  ripple  is  used  the  process  of 
reversal  is  unnecessary. 

A  velvet- like  black  coating  of  finely 

divided  platinum  may  also  be  obtained  by  chemical  immersion  techniques,  i.e.. 
by  using  a  metal  more  active  than  platinum  to  replace  the  platinum  from 
solution.  Black  coatings  were  obtained  on  nickel  in  this  manner. 

13.2.2.10  Carbon  Black  on  Aluminum 

A  deep  black,  uniform  coating  can  be  obtained 
on  aluminum  by  depositing  carbon  from  the  incomplete  combustion  of  natural 
gas.  This  system  is  of  interest  because  of  its  ease  of  application  and  the 
ease  with  which  its  reflectivity  can  be  varied  over  wide  limits.  On  smooth 
surfaces,  its  adhesion  is  not  good,  but  the  adhesion  is  much  improved  on  a 
microscopically  rough  surface. 

Carbon  has  been  deposited  on  dull  and  bright 
platinum,  copper,  and  aluminum.  The  most  successful  of  the  carbon  black 
systems  which  our  laboratory  has  investigated  is  the  carbon  black  on  alum¬ 
inum.  The  procedure  for  obtaining  this  surface  is  as  follows: 

1.  Preparation  of  aluminum  base: 

The  aluminum  sheet  which  was  used  was  obtained  from 
the  Aluminum  Company  of  America  and  is  their  Standard 
Bright  Finish,  Grade  1100,  H-14.  The  aluminum  sheet 
did  not  require  any  preliminary  cleaning  treatment. 

2.  Preparation  of  surfaces: 

The  aluminum  sheet  test  sample  (2"  x  2")  was  held 
for  about  five  minutes  in  the  tip  of  an  8”-  9”  luminous 
flame  from  a  Fisher  burner  which  had  its  air  supply 
completely  sealed  off.  Only  enough  carbon  is 
deposited  to  give  the  surface  a  rust-colored  cast 
when  viewed  orthogonal  to  the  plate.  The  surface 
has  a  deep  black  cast  when  viewed  nearly  parallel  to 
the  plate. 
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The  thickness  of  the  aluminum  oxide  film 


which  forms  on  the  aluminum  surface  will  be  a  function  of  the  plate  temper¬ 
ature,  and  monitoring  the  temperature  of  the  plate  with  a  thermocouple  could 
serve  as  an  additional  experimental  control.  This  system  has  not  been  fully 
evaluated  because  of  difficulties  in  measuring  the  carbon  film  thickness. 

13.2.2.11  Electrically  Conducting  Metal  Oxide  Coatings 

on  Glass  (NESA  Glass) 

Many  procedures  exist  in  the  literature  for 
preparing  these  conductive  glasses  (Refs.  26  and  27),  which  generally 
exhibit  low  reflection  in  the  visible  and  high  reflection  in  the  infrared. 

They  could  therefore  act  as  plates  used  to  retain  incident  solar  radiation 
within  a  cavity.  Corning  Glass  Works  provides  commercially  an  infrared 
reflecting  glass  which  reflects,  according  to  their  literature,  from  6-8 
percent  in  the  visible,  30  percent  at  2.2  microns,  and  approaches  70  percent 
reflection  as  an  asymptote  at  5  microns.,  Samples  have  been  received  and 
evaluated  in  this  laboratory. 

The  original  EOS  attempts  at  producing 
electrically  conducting  glass  coatings  were  made  by  spraying  a  stannic 
chloride- isopropyl  alcohol  solution  from  an  atomizer  onto  a  heated  glass 
surface.  Both  soft  glass  and  fused  quarts  were  used  and  an  optimum  tem¬ 
perature  was  found  at  600  -  700°C.  The  minimum  resistance  for  these  coatings 
was  200  ohms  per  square.  For  optimum  spectral  characteristics  the  resistance 
should  be  about  5  -  20  ohms  per.  square.  In  another  method,  SnCl^ 

\>as  put  into  a  heated  bulge  in  a  long  Vycor  blowpipe.  Over  this  was  passed 
air  saturated  with  water  vapor  while  the  blowpipe  was  held  in  a  horizontal 
position.  From  the  tip  of  the  tube  the  vapor  impinged  upon  the  glass. 

Optimum  conditions  were  found  for  the  temperature  of  the  SnCl^  *  SH^O, 
the  temperature  of  the  glass,  and  for  the  air -water  ratio,  which  was  varied 
by  changing  the  temperature  of  the  water  used  in  saturating  the  air.  The 
lowest  'resistance  obtained  in  either  of  these  methods  was  200  ohms  per  square. 

The  hydrolysis  of  stannic  chloride  is  very  slow 
at  room  temperature  and  occurs  mostly  on  the  hot  surface  of  the  glass. 
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However,  there  is  a  large  amount  of  hydrolysis  which  occurs  in  the  heated 
liquid  SnCl^  •  5H^0  in  the  blow  pipe.  In  order  to  prevent  this,  calcium 

chloride  was  mixed  with  the  SnCl^  •  5^0  as  a  ^ry^nS  agent.  In  this 

manner  a  superior  coating  with  a  resistance  of  150  ohms  per  square  was 
obtained.  Pure  anhydrous  stannic  chloride  was  also  tried  but  was  given 
up  since  a  resistance  of  410  ohms  per  square  was  the  lowest  obtainable 
for  very  short  duration  spray  times. 

Solution  No.  2  from  Patent  No.  2,564,706 
(Ref.  13)  by  John  M.  Mochel  (assigned  to  the  Corning  Glass  Works)  was  then 
tried.  This  proved  very  successful  since  it  was  possible,  by  this  method, 
to  get  resistances  from  less  than  five  to  about  twenty  ohms  per  square. 

The  solution  contained  50  ml.  of  water,  10  ml.  of  hydrochloric  acid, 

100  gm.  of  stannic  chloride  pentahydrate  (SnC14  •  5H^O),  and  0.5  gm.  of 

antimony  trichloride  (SbCl^),  It  was  put  into  a  DeVilbis  No.  15  aspirator 
through  which  air  was  passed  from  a  compressor.  A  pressure  of  about  45  cm. 
of  mercury  was  found  necessary  to  produce  the  proper  drop  size  for  a  uni¬ 
form  coating.  For  these  purposes  the  optimum  temperature  was  found  to  be 
between  600  and  650°C  for  the  surface  of  the  glass.  The  glass  plates  were 
laid  flat  on  the  floor  of  a  muffle  furnace  and  sprayed  for  about  one 
minute  from  a  distance  of  about  one  foot,  with  the  spray  directed  onto  the 
surface  at  a  slight  angle. 

The  cleaning  of  the  glass  plates  is  an 
important  factor.  The  stannic  oxide  coating  of  glass  plates  which  were 
uncleaned  had  a  resistance  of  1500  ohms  per  square;  whereas  those  cleaned 
with  sulfuric  ac id -di chromate  cleaning  solution  had  a  resistance  of  five 
to  twenty  ohms  per  square. 

The  relationship  between  surface  resistance 
and  spectral  properties  is  discussed  in  Section  13.2.3. 

13.2.3  Spectral  Evaluation  of  Surfaces 

Each  selective  surface  sample  produced  in  the  EOS  labor¬ 
atory  was  tested  in  the  Perkin-Elmer  13U  spectrophotometer  in  order  to 
determine  reflectance  or  transmittance  as  a  function  of  wavelength.  This 
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instrument  is  described  in  Section  13.5.3. 

The  procedures  used  to  record  the  reflectance  of  both 
the  black  coated  surfaces  and  the  NESA  glass  were  essentially  the  same. 
These  measurements  are  not  absolute,  but  are  based  on  the  known  reflectance 
characteristic  of  an  aluminum  mirror.  Therefore,  after  standardizing  the 
reflectance  scale  in  terms  of  the  aluminum  characteristic,  the  aluminum 
was  replaced  by  the  actual  samples  and  their  reflectances  determined  over 
the  same  range. 

Gold  smoke  transmittance  curves  were  obtained  directly 
in  the  double  beam  mode  of  operation  of  the  spectrophotometer.  This  allows 
an  absolute  transmittance  curve  to  be  recorded  directly  on  the  chart.  A 
study  was  conducted  to  determine  what  effect  changing  the  direction  through 
which  the  light  was  passed  through  the  NESA  glass  had  on  the  transmittance 
in  the  visible  region  of  the  spectrum.  Figure  36  compares  the  transmittance 
measurements  when  the  coated  side  was  facing  the  light  source  and  when  the 
coated  side  is  away  from  the  light  source.  On  the  basis  of  this  study  it 
can  be  concluded  that  the  transmission  of  NESA  glass  is  unaffected  by  the 
direction  through  which  light  is  passed. 

One  generalization  which  can  tentatively  be  made  for  the 
opaque  spectrally  selective  surfaces  produced  at  Electro-Optical  Systems 
is  that  the  solar  absorption  and  infrared  emission  both  increase  as  the 
thickness  of  the  coating  increases  and  that  for  certain  surfaces  the 
inflection  point  shifts  toward  the  longer  wavelengths.  The  absorption 
and  emission  can  thus  be  controlled  to  some  extent,  but  not  independently., 
Furthermore,  the  quantitative  functional  dependence  of  a  and  €  on  coating 
thickness  is  not  now  known. 

This  section  contains  brief  descriptions  of  the  spectral 
characteristics  of  each  surface  produced.  In  some  cases,  conclusions  are 
drawn  regarding  the  reasons  for  the  results  obtained. 


92 


1. 


Gold  Smoke  on  Glass 


The  results  of  the  gold  smoke  investigations  are 
given  in  Fig.  34.  It  is  clear  from  these  curves  that  widely  varying 
spectral  characteristics  may  result  from  small  changes  in  deposition 
procedure.  Samples  with  good  visible  absorption  also  tended  to  show 
relatively  poorer  infrared  transmission. 

The  better  samples  developed  during  this  program 
showed  spectral  characteristics  generally  in  accordance  with  those  obtained 
by  other  workers  (Refs.  17,18,  and  19).  However,  considerably  improved 
coatings  may  be  possible  by  means  of  techniques  currently  being  investigated 
at  EOS  on  other  programs. 

2 .  Sprayed  Cupric  Oxide  on  Aluminum 

The  results  obtained  from  these  surfaces  (illustrated 
in  Fig.  37)  agreed  well  with  those  obtained  by  Unger  (Ref.  13  and  Fig.  38). 

The  best  of  the  Unger  surfaces  was  used  in  the  heat  equilibrium  tests.  The 
surface  was  examined  in  two  areas.  The  reflectance  data  for  these  two  areas 
are  compared  in  Fig.  37  and  indicate  the  degree  of  uniformity  of  the  surface. 
Inasmuch  as  Unger  carried  out  a  careful  and  thorough  research  program  on  these 
coatings,  it  is  unlikely  that  they  will  be  greatly  improved  in  the  future. 
However,  because  of  the  research  already  expended,  they  are  among  the  best 
spectrally  selective  surfaces  yet  produced  in  regard  to  both  their  spectral 
and  physical  properties. 

The  only  additional  comment  which  our  laboratory  has 
to  make  regarding  the  fabrication  of  Unger  surfaces  is  that,  in  general, 
the  spray  rate,  the  spray  time,  and  the  method  of  spraying  must  be  adjusted 
so  that  the  baked  coat  is  no  darker  than  a  very  light  gray.  Evidently  the 
film  must  be  thin  enough  to  be  transparent  in  the  infrared  and  allow  the  low 
emissivity  characteristics  of  the  aluminum  to  be  the  controlling  factor.  A 
reflectance  versus  wavelength  curve  (Fig.  39)  for  uncoated  99,99  percent  pure 
aluminum  is  included  so  that  a  comparison  can  be  made  between  the  coated  and 
uncoated  surfaces.  Aluminum  Sample  No.  2  was  used  in  the  preparation  of  this 
Unger  surface. 
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A  great  deal  of  effort  was  spent  in  the  selection  of 
the  proper  aluminum  substrate.  The  following  types  of  aluminum  surfaces 
were  investigated: 

Sample  No.  1  99.99  percent  aluminum,  0.662"  thick,  1/2-H 

Sample  No.  2  99.99  percent  aluminum,  0.062”  thick,  1/2- H 

Sample  No.  3  99.99  percent  aluminum,  0.062"  thick,  1/2-H 

Sample  No.  4  Reflector  grade  aluminum,  0.040"  thick 

Sample  No.  1  and  Sample  No.  2  were  from  the  same  sheet  of  aluminum,  only 
Sample  No.  1  was  taken  from  a  portion  of  the  seet  which  appeared  to  have  the 
least  visible  reflectivity.  Sample  No.  2  was  taken  from  a  portion  of  the 
sheet  which  appeared  to  have  the  greatest  visible  reflectivity.  Sample  No,  3 
was  buffed  on  a  buffing  wheel.  Sample  No.  4  is  what  the  Aluminum  Company  of 
America  calls  their  reflector  grade  aluminum.  Figure  39  gives  a  comparison 
of  the  reflectance  between  these  four  samples.  Sample  No.  2  was  chosen  to 
be  used  in  the  preparation  of  the  Unger  surface.  Another  sample  of  aluminum 
sheet  which  was  evaluated  for  use  as  an  Unger  substrate  was  99.99  percent 
pure  aluminum,  0.01  inch  thick,  0-H  temper.  This  aluminum  has  a  uniformly 
bright  surface  to  the  eye.  The  0.01  inch  thick  aluminum  was  not  spectrally 
evaluated  with  the  other  samples  because  it  was  too  thin  to  be  used  £n  the 
heat  equilibrium  test  unit.  The  99.99  percent  aluminum,  0.062  inch  thick, 
is  1/2-H  temper  and  has  a  dull  coating  in  places  (probably  a  thick  oxide  film) 
which  may  have  been  introduced  during  the  tempering  process. 

Aluminum  Sample  No.  2  was  preferred  to  aluminum  Sample  No.  3 
because  it  was  felt  that  the  compound  used  in  buffing  might  not  be  entirely 
removed  by  the  organic  solvents  and  that  this  might  introduce  another  experimental 
variable  which  would  require  control. 

3 .  Cupric  Oxide  on  Copper  by  the  "Ebonol"  Process 

As  shown  in  Fig.  40,  this  system  has  0  percent  reflectance 
from  0.3  -  6  microns  and  becomes  reflecting  in  the  range  6-9  microns.  No 
conclusions  can  be  made  as  to  the  difference  in  reflectance  between  Samples  0, 

1,  2,  and  3,  all  of  which  were  prepared  by  the  same  procedure.  The  reflectance 
is  less  than  50  percent;  therefore,  without  further  investigation  it  does  not 
appear  promising  for  the  present  application.  Further  Investigations  should 
Include  forming  a  coat  as  thin  as  possible  on  a  highly  polished  substrate. 
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Polished  copper  was  spectrally  evaluated,  and  found  to 
be  as  reflective  in  the  infrared  as  aluminum.  Etched  copper  has  a  specular- 
reflectance  much  lower  than  that  of  polished  copper;  it  may  or  may  not  have 
a  lower  diffuse  reflectance. 

Although  it  was  desirable  to  prepare  and  evaluate  thin 
cupric  oxide  coatings  on  polished  copper,  this  was  not  done  since  the  copper 
substrate  had  to  be  etched  prior  to  oxidation  and,  therefore,  the  coated  surface 
could  not  be  evaluated  by  specular  reflectance. 

Cupric  oxide  coatings  using  Ebonol  have  been  prepared  and 
evaluated  by  other  workers  in  this  field  (Ref.  28)  Gier  and  Dunkle  claim  to 
have  produced  surfaces  with  five  percent  reflectance  in  the  visible  and  80  -  85 
percent  reflectance  in  the  infrared.  The  inflection  point,  however,  is  at  only 
one  micron  which  would  give  a  trapping  efficiency  of  only  70  percent  for  an 
equilibrium  temperature  of  400°C. 

4.  Black  Chromate  on  Magnesium 

This  coating  was  powdery  and  non-adherent  so  spectral 
measurements  were  not  conducted. 

5 ♦  Black  Nickel-Iron  Oxide  on  Steel  by  the  "Nickel 

Pent r a te"  Process  (see  Fig.  41) 

These  surfaces  had  0  percent  reflectance  from  0.3  -  2.25 
microns.  The  surface  begins  to  reflect  at  2.25  microns.  At  8.5  microns  the 
reflectance  is  17.5  percent.  The  spectral  characteristics  of  this  system  are 
not  satisfactory  for  this  application,  so  no  further  work  was  done  on  this 
surface. 

6.  Black  Nickel  on  Aluminum 

These  surfaces  had  0  percent  reflectance  from  0.9  microns 
to  7.5  microns.  The  spectral  characteristics  of  this  system  are  also  not 
satisfactory  for  our  purpose  so  no  further  work  was  done  on  this  system. 

7 .  Black  Nickel  on  Copper 

As  shown  in  Fig.  42,  this  surface  is  one  of  the  most 
promising  for  the  present  application.  It  has  less  than  10  percent  reflectance 
in  the  visible  and  is  greater  than  95  percent  reflecting  in  the  infrared  with 
a  very  sharp  cutoff  between  1.0  and  2.5  microns.  The  best  of  the  black  nickel 
surfaces  was  used  in  the  heat  equilibrium  tests.  The  surface  was  examined  in 
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two  areas,  diagonally  across  the  plate  from  each  other.  The  reflectance 
data  for  these  two  areas  are  compared  in  Fig.  42  and  indicate  the  degree 
of  uniformity  of  the  surface. 

In  the  course  of  studying  the  black  nickel  coating 
there  was  noted  a  considerable  difference  in  optical  characteristics 
between  samples.  (See  Fig.  43).  Three  variables  in  the  procedure  (time 
of  coating ,  current  density,  and  bath  temperature)  were  controlled  so  as 
to  determine  their  effect  on  the  optical  characteristics.  Fourteen  samples 
were  prepared  while  the  experimental  variables  were  changed  within  the 
recommended  limits: 

Time  of  coating  1-3  minutes 

2 

Current  density  1-2  amps/ft 

Bath  temperature  27  -  38°C 

It  was  observed  that  at  lower  temperatures  (less  than  36°C),  and  as  the 
bath  begins  to  become  depleted,  the  casting  received  a  definitely  yellow 
cast.  It  was  also  determined  that  voltage  control  was  extremely  important, 

1.5  volts  being  the  optimum.  Both  current  and  voltage  should  be  metered 
during  the  electrolytic  process.  Undetermined  conditions,  which  cause  the 
surface  to  appear  blue-black,  also  seem  to  produce  surfaces  with  improved 
spectral  characteristics  for  this  application.  Further  study  is  required 
in  order  to  determine  the  exact  nature  of  the  conditions  which  produce 
variations  in  the  spectral  response  of  these  surfaces. 

8.  Black  Nickel  on  Nickel  and  Nickel  Plated  Copper 

This  surface  is  also  very  promising  as  shown  in  Fig.  44. 
Toward  the  close  of  the  experimental  work,  it  was  concluded  that  in  the 
evaluation  of  the  spectral  characteristics  of  surfaces,  it  would  be  helpful 
to  determine  the  reflectance  of  the  uncoated  plate  before  the  plate  receives 
surface  treatment.  The  inflection  point  on  the  curve  of  reflectance  versus 
wavelength  for  several  of  these  black  surfaces  can  be  shifted  toward  the 
visible  merely  by  polishing  the  substrate.  A  comparison  of  the  results  of 
the  black  nickel  on  nickel  system  and  the  black  nickel  on  nickel-plated 
copper  illustrates  this  point.  See  Fig.  44.  The  black  nickel  was  deposited 
on  nickel  polished  to  a  near  mirror-bright  finish  and  bn  nickel-plated  copper, 
which  was,  prior  to  plating,  polished  to  a  mirror-bright  finish.  The  effect 
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of  the  brighter  base  on  the  shift  of  the  inflection  point  toward  the  visible 
is  evidenced  here . 

9 .  Cupric  Oxide  on  Copper  by  the  "Dura-Black"  Process 
This  surface  also  showed  promise  for  solar  energy 
collection.  As  illustrated  in  Fig.  31,  the  system  is  0  percent  reflecting 
in  the  visible,  begins  to  reflect  at  1  micron  and  becomes  100  percent 
reflecting  in  the  infrared. 

10 .  Platinum  Black  on  Platinum 

To  our  knowledge,  platinum  black  coatings  on  platinum 
have  never  before  been  evaluated  as  potential  spectrally  selective  surfaces. 

Fig.  45  shows  the  reflectance  of  platinum  black  on  a  platinum  surface  as  a 
function  of  wavelength.  This  surface  is  mechanically  unstable;  however,  it 
nas  tne  unique  adv  antage  of  chemical  inertness  and  thermal  stability,  which 
would  be  very  useful  in  a  system  in  which  these  qualities  would  be  very 
important.  This  surface  has  probably  the  best  spectral  characteristics  for 
the  purpose  of  trapping  solar  energy  of  any  of  the  surfaces  prepared  and 
investigated  by  our  laboratory.  Tests  have  been  made  which  indicate  that  the 
platinum  black  film  thickness  controls,  to  some  extent,  the  inflection  point, 
the  degree  of  slope,  and  the  infrared  emission.  With  Unger* s  surface,  in 
order  to  gain  low  emissivity  in  the  infrared,  absorptivity  in  the  visible  has 
to  be  sacrificed.  In  the  case  of  platinum  black  on  platinum  there  is  no  need 
to  sacrifice  in  the  visible  to  gain  in  the  infrared. 

Platinum  black  has  been  electrochemically  deposited  on 
platinum- plated  copper;  however,  the  substrate  metal  was  active  and  contaminated 
the  bath  thus  ruining  the  coating.  In  the  future  there  should  be  an  invest¬ 
igation  of  platinum  black  on  a  thick  electroplated  platinum  plate. 

The  mechanical  stability  of  the  platinum  black  coating 
may  be  increased  by  converting  the  platinum  black  to  a  more  adherent  platinum 
coating,  .namely  platinum  gray.  This  can  be  accomplished  by  heating  the 
platinum  black  coated  surface  to  a  dull  redness.  This  conversion  may  or  may 
not  alter  the  optical  properties  of  this  system.  The  adhesion  of  platinum 
black  to  platinum  may  also  be  increased  by  microscopically  roughening  the 
platinum  substrate. 
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Platinum  black  has  been  chemically  plated  on  nickel, 
aluminum,  and  copper.  Whereas  platinum  black  on  platinum  is  not  very 
adherent,  platinum  black  on  nickel  and  copper  is  very  adherent.  Platinum 
black  on  nickel  has  been  spectrally  evaluated  and  does  not  appear  promising. 
However,  the  nickel  substrate,  which  has  not  been  spectrally  evaluated,  appears 
dull  to  the  eye  and  this  may  account  for  the  poor  spectral  results..  Further 
effort  in  this  phase  of  the  work  is  advised.  Platinum  black  on  copper  appears 
very  promising.  It  is  too  early  to  judge  the  promise  of  platinum  black  on 
aluminum. 

11 .  Carbon  Black  on  Aluminum 

To  our  knowledge,  carbon  black  coatings  on  highly 
reflective  metals  have  never  before  been  evaluated  as  potential  spectrally 
selective  surfaces.  Carbon  black,  though  mechanically  unstable,  offers 
exceptional  absorptance  in  the  visible  with  a  good  control  of  emittance 
in  the  infrared.  See  Fig.  46.  It'  has  been  established,  as  in  the  case  of 
platinum  black,  that  the  film  thickness,  to  some  degree,  controls  the  inflection 
point,  the  degree  of  slope,  and  the  infrared  emission.  Although  film  thickness 
could  not  be  quantitatively  controlled,  it  could  be  qualitatively  determined  by 
visual  observation  and  by  noting  the  coating  time.  Surfaces  with  increasingly 
thick  deposits  had  increasing  absorption  in  the  visible,  a  decrease  in  slope, 
an  inflection  point  shift  toward  the  infrared,  and  an  increase  in  infrared 
emission.  It  is  strongly  felt  that  the  spectral  characteristics  of  this 
surface  can  be  improved  by  depositing  extremely  thin  carbon  films. 

12 .  Metal  Oxide  Coat  infs  on  Glass  (NESA  Glass) 

The  characteristic  desired  with  metal  oxide  coated  glass 
is  high  solar  transmittance,  and  high  infrared  reflectance.  Therefore,  trans¬ 
mittance  as  well  as  reflectance  is  determined  for  the  glass.  For  the  purposes 
of  this  program,  the  optical  properties  of  thin  pyrex  glass  are  satisfactory. 

The  metal  oxide  conductive  coatings  have  also  been  produced  by  Electro-Optical 
Systems  on  fused  silica  test  pieces  at  700°C  with  identical  results.  The 
advantage  to  be  gained  through  the  use  of  fused  silica  is  a  slightly  higher 
transmittance  for  the  same  thickness  of  glass.  Spectrally,  the  transmittance 
of  both  types  of  glass  is  flat  from  0.3  to  2.3  microns. 
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The  characteristics  of  six  samples  of  NESA  glass  which 
were  the  best  prepared  at  our  laboratory  and  one  prepared  at  Corning  Glass 
Works  are  shown  in  Figs.  32 ,  33 ,  47,  48 ,  and  49.  In  Fig.  47,  it  is  seen 
that  NESA  Sample  No.  29  approaches  both  in  electrical  resistance  and  reflect¬ 
ance  very  closely  to  that  of  Corning  Glass.  NESA  Samples  No.  27,  28,  and  29 
were  prepared  on  1/16”  soft  glass;  NESA  Samples  P-1  and  P-2  (Figs.  48  and  49) 
were  prepared  on  1/8”  Pyrex  glass,  and  NESA  Sample  No.  P-6  (Figs.  32  and  33) 
was  prepared  on  1/16"  Pyrex  glass.  NESA  Sample  No.  P-6  was  selected  for  use 
in  the  heat  equilibrium  tests.  Aside  from  the  difference  in  glass  base,  the 
chemical  formulation  of  the  solution  used  in  spraying  and  the  method  of 
preparation  is  the  same  for  each  of  the  above  NESA  samples. 

One  generalization" which  can  tentatively  be  made  for  the 
transparent  spectrally  selective  surfaces  produced  at  EOS  is  in  regard  to  the 
relationship  between  surface  resistance  and  spectral  properties.  Coatings 
prepared  with  stannic  chloride  hydrate  and  pure  anhydrous  stannic  chloride 
had  resistances  ranging  from  150  to  1,300,000  ohms  per  square;  while  coat¬ 
ings  prepared  with  solution  No.  2  from  Patent  No.  2,564,706  had  resistances 
ranging  from  five  to  twenty  ohms  per  square.  The  higher  resistance  coatings 
had  spectral  characteristics  which  were  unsuitable  for  use  in  a  solar  heat 
trap,  whereas,  the  lower  resistance  coatings  have  spectral  characteristics 
which  are  very  satisfactory  for  this  purpose. 

The  presence  of  antimony  trichloride  is  necessary  to 
produce  surfaces  of  low  resistance.  Chemical  compounds  other  than  antimony 
trichloride  have  been  used  to  produce  surfaces  of  low  electrical  resistance. 
Patent  2,694,649  (1954)  indicates  that  an  In^O^  coating  on  glass  can  give  as  low 
a  resistance  as  the  SnO^,  Sb^O^  coat*n8s  Prepared  from  the  Corning  patent.  It 
would  be  interesting  to  compare  the  spectral  properties. 

Patent  2,756,165  (1956)  describes  a  process  in  which 
transparent  plastic  is  coated  with  a  material  having  high  transmission  and 
low  reflection  in  the  visible.  It  would  be  interesting  to  determine  if  the 
reflectance  in  the  infrared  is  high  enough  to  b&  of  interest. 

In  order  to  exactly  determine  the  optimum  resistance  for 
metal  oxide  coatings,  a  comparison  must  be  made  between  energy  lost  from  low 
transmittance  in  the  visible  and  energy  gained  from  high  reflectance  in  the 
infrared;  i.e.,  both  reflectance  and  transmittance  must  be  optimized  as 
functions  of  resistance. 
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The  reflectance  of  the  visible  part  of  the  spectrum 
seems  to  be  connected  in  no  simple  fashion  with  the  resistance  of  those 
coatings  which  have  resistances  in  the  range  of  five  to  twenty  ohms  per 
square.  Figure  50  is  a  plot  of  reflectance  versus  wavelengths  for  NESA 
glasses  with  resistances  ranging  from  4.9  to  20.0  ohms  per  square. 

The  reflectance  in  the  infrared,  however,  shows  a 
definite  correlation  with  the  resistance  of  the  coated  surface.  Figure 
47  shows  reflectance  characteristics  for  NESA  glasses  with  resistances 
from  6.6  to  20.0  ohms.  The  surfaces  with  electrical  resistances  of  6.6 
and  11.0  ohms  pel  square  have  almost  10  percent  higher  reflectance  than 
do  surfaces  with  resistances  of  14.7  and  20.0  ohms  per  square.  This 
correlation  of  electrical  resistance  and  infrared  reflection  is  also  borne 
out  by  another  set  of  stannic  oxide  glasses.  P-1  and  P-2  (see  Fig.  48), 
which  were  prepared  here  in  our  laboratories'.  NESA  Sample  No.  P-1  which 
has  a  resistance  of  16.0  ohms  per  square  has  a  15  percent  lower  reflectance 
in  the  infrared  than  does  NESA  Sample  No;  P-2  with  a  resistance  of  only 
3.7  ohms.  Based  on  limited  experimental  data  it  can  be  concluded  that  the 
lower  the  electrical  resistance  the  higher  the  reflectance  in  the  infrared 
portion  of  the  spectrum. 

The  transmittance  in  the  visible  part  of  the  spectrum 
shows  a  correlation  with  electrical  resistance  of  the  surface  coating,  see 
Fig.  49.  P-1  with  a  resistance  of  16  ohms  per  square  has  a  transmittance 

about  30  percent  higher  than  P-2  with  a  resistance  of  2.7  ohms  per  square. 

In  Fig.  36,  a  comparison  is  made  between  the  trans¬ 
mittance  of  uncoated  Pyrex  glass  and  Pyrex  glass  coated  with  stannic  oxide. 
The  stannic  oxide  coating  decreases  the  transmittance  cf  the  pyrex  glass  in 
the  visible  region  of  the  solar  spectrum.  This  loss  may  or  may  not  be  offset 
by  the  increased  reflectance  of  the  coated  glass  in  the  infrared  region, 
where  the  J lack  surface  is  emitting,  This  comparison  of  the  effectiveness  of 
NESA  versus  uncoated  glass  is  left  for  future  analytical  and  experimental 
s  tudy . 
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13.2.4  Conclusions  Regarding  Suitability  of  Surfaces 

The  suitability  of  the  various  spectrally  selective 
surfaces  for  solar  energy  collection  depends  on  their  spectral  character¬ 
istics  ,  uniformity,  ease  of  fabrication,  and  cost.  A  comparison  of 
properties  for  the  best  opaque  spectrally  selective  surfaces  prepared  at 
Electro-Optical  Systems  is  given  in  Table  IV..  The  description  of  each 
property  is  based  only  on  a  qualitative  comparison  between  the  systems 
listed  and  not  upon  any  absolute  standard.  The  solar  reflectivity  gives 
the  percent  reflectance  in  the  visible  portion  of  the  spectrum.  The 
infrared  emissivity  gives  the  percent  emittanee  from  7  -  15  microns. 

Of  the  systems  which  have  been  evaluated,  the  best, 
based  on  over-all  considerations,  seem  to  be  the  black  nickel  on  nickel 
or  copper  and  the  cupric  oxide  on  aluminum  by  the  Unger  process. 


13.3  Radiation  Traps 

A  radiation  trap  may  be  defined  as  a  device  which  exhibits  a 
high  absorptance  for  solar  energy  while  displaying  the  characteristics 
of  a  poor  emitter  in  the  long  wave  length  region.  There  are  two  basic 
schemes  which  may  be  used  to  realize  the  trapping  effect;  one  differen¬ 
tiates  between  the  wavelengths  of  radiation,  while  the  other  distinguishes 
between  the  wavefronts  of  incident  and  reradiated  energy  by  geometric  means. 

A  surface  having  the  former  characteristic  is  termed  spectrally  selective. 

The  spectrally  selective  trap  consists  of  a  specially  treated 
surface  which  absorbs  solar  radiation  but  is  a  poor  emitter  for  the  long 
wavelength  thermal  radiation.  Several  methods  for  obtaining  such  surfaces 
are  described  elsewhere.  (See  Section  13.2).  As  an  example,  chemical 
and  electrolytic  deposition  of  thin  black  coatings  on  polished  aluminum 
plates  yield  such  surfaces.  Another  approach  is  the  treatment  of  cover 
glasses  so  that  they  will  transmit  solar  radiation  to  a  black  absorber  plate 
beneath  but  will  reflect  back  the  thermal  radiation  from  absorber. 

Wavefront  discrimination  is  accomplished  by  geometric  means.  A 
familiar  example  of  such  a  trap  is  the  use  of  an  insulated  cavity  receiver 
in  conjunction  with  a  concentrator  (lens  or  mirror).  The  parallel  solar 
radiation  is  focussed  through  a  small  hole  in  the  cavity  receiver  and  almost 
totally  absorbed  on  the  walls  of  the  cavity.  Only  a  very  small  amount  of  the 
diffuse  radiation  from  the  inside  walls  of  the  cavity  can  escape  from  the 
small  hole., 

The  analytical  and  experimental  results  from  the  collector  program, 
as  well  as  the  over-all  system  requirements,  have  indicated  a  heat  collector- 
exchanger  design  which  may  satisfy  the  requirements  of  the  solar  propulsion 
system.  This  design  (see  Section  11.3)  consists  of  a  series  of  airfoil 
shaped  struts  placed  parallel  to  the  airflow  in  the  heat  exchanger.  Incoming 
radiation  is  focussed  into  the  strut  through  a  hole  in  the  top  surface  by 
means  of  a  Fresnel  lens.  After  passing  through  the  hole  the  energy  is  scattered 
within  the  strut  and  absorbed  on  the  inside  surface.  Energy  losses  in  this 
system  can  only  occur  as  reradiation  from  the  surface  of  the  hole,  or  from 
losses  in  the  focussing  system.  These  losses  are  discussed  in  the  first 
subsection.  The  second  subsection  discusses  the  performance  of  traps  made 
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and  tested  in  the  Electro-Optical  Systems,  Inc.  laboratory. 

13.3.1  Energy  Losses  from  Geometric  Radiation  Traps 

In  evaluating  the  energy  losses  of  geometric  heat  traps 
it  is  necessary  to  consider  the  transmission  efficiency  of  the  lens-entrance 
hole  combination  and  the  heat  loss  by  radiation  from  the  hole.  These  factors 
are  clearly  interdependent  since  enlarging  the  hole  increases  the  transmission 
efficiency  but  also  increases  the  heat  loss  from  the  hole.  The  following 
will  refer  specifically  to  the  airfoil-shaped  geometric  radiation  trap- 
However,  most  of  the  comments  are  generally  applicable. 

Suppose  the  strut  has  a  total  surface  area,  Ag>  emissivity 
€g,  and  is  at  a  uniform  temperature,  T.  If  the  strut  is  radiating  to  an 
environment  whose  equivalent  temperature  is  sufficiently  low,  the  heat  loss 
by  radiation  from  the  strut  surface,  Qg,  may  be  written: 

Qs  =  Ag  €s  CT  T4  (263) 


If  the  strut  entrance  hole  has  a  diameter  d,  and  is 

n 

radiating  like  a  black  body  at  temperature  T,  the  heat  loss  from  the 
hole  to  a  low  temperature  environment  is  given  by: 


a  T4 


(264) 


The  diameter  of  the  entrance  hole  will  probably  be  made  several  times 
larger  than  the  image  diameter  of  the  sun  in  order  to  account  for 
aberrations,  lens  defects  and  misalignments.  Let  |  be  this  factor. 

If  f  is  the  focal  length  of  the  lens,  then  the  hole  diameter  may  be  written 
approximately  as: 


because  the  sun  (and,  thus,  the  image)  subtends  an  angle  of  about  1/100 
radian  as  seen  from  the  lens.  Substituting  this  into  the  last  equation 
and  taking  the  ratio  Qh/Qg  yields:  \ 
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(266) 


2.2 

± JJL 


4  x  10  A  € 
s  s 


If  the  lens  is  circular  with  a  diameter  and  a  projected  area  A^, 
the  above  equation  may  be  written 


S 

% 


(267) 


Either  of  the  last  two  equations  may  be  evaluated  in  order  to  predict 
the  importance  of  the  heat  losses  through  the  entrance  hole.  If  the 
design  parameters  are  taken  to  be  as  follows: 


f 


1 

2 


1 

2 


(equal  projected  area) 


and  an  optimistic  value  of  0.02  is  chosen  for  eg  (which  tends  to  increase 
Qh/Qs);  then  even  for  £  =  10,  the  heat  loss  from  the  entrance  hole  is  only 
6  percent  of  the  losses  from  the  airfoil  surface.  The  assumption  made  that 
the  heat  loss  from  the  entrance  hole  is  negligible  appears  to  be  justified. 

In  evaluating  the  transmission  efficiency  of  the  lens-entrance 
hole  combination  it  is  necessary  to  know  the  distribution  of  flux  at  the  focus 
of  the  lens.  Once  this  distribution  is  known,  the  blocking  effect  of  small 
entrance  holes  may  be  evaluated.  Clearly,  if  the  entrance  hole  is  made  large 
enough,  the  transmission  efficiency  of  the  combination  is  simply  that  of  the 
lens  itself.  Thin  sheet  plastic  Fresnel  lenses  have  been  made  which  exhibit 
a  transmission  of  75  to  85  percent. 

Depending  upon  the  emissivity  of  the  strut  surface,  entrance 
holes  whose  diameters  are  10  to  20  times  larger  than  the  solar  image  may  be 
used  without  introducing  serious  radiation  losses.  (Equations  266  ana  267). 
Hence  it  is  reasonable  to  assume  transmission  efficiencies  of  the  lens-entrance 
hole  combination  of  about  70  to  80  percent.  The  exact  values  will,  of  course, 
depend  upon  the  final  design  of  the  lens  including  the  mechanical  perfection 
achieved. 
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In  order  to  reduce  the  heat  losses  from  entrance  holes, 

It  is  often  suggested  that  NESA  glass  covers  be  used  over  the  holes.  Light 
passes  through  the  NESA  glass  and  is  absorbed  by  the  "black"  surface  which 
reaches  a  temperature  T.  Part  of  the  energy  which  is  reradiated  from  the 
inside  of  the  trap  structure  impinges  upon,  and  is  mostly  reflected  by, 
the  NESA  glass.  That  component  which  is  not  reflected  is  absorbed  and 
serves,  along  with  the  small  amount  of  incoming  light  energy  which  is 
absorbed,  to  raise  the  temperature  of  the  glass.  The  NESA  glass,  in  turn, 
radiates,  to  both  the  inside  of  the  trap  and  to  the  "sky".  This  upward 
loss  of  energy  is  less  than  that  which  would  be  lost  through  the  hole  itself 
since  the  NESA  glass  is  at  a  lower  temperature  than  the  strut.  However,  a 
new  energy  loss  has  been  introduced;  approximately  30  percent  of  the  light 
energy  striking  the  NESA  glass  is  reflected.  (Various  coatings  can  be  used 
to  reduce  reflection  losses  slightly).  Thus,  it  is  only  for  very  large 
holes,  where  the  radiation  losses  would  normally  be  somewhat  more  than 
30  percent  of  the  input  energy,  that  using  the  NESA  glass  in  this  manner 
pays  off.  The  crossover  point  can  be  calculated  as  follows: 

Since  the  area  of  the  NESA  window  is  small  compared  with 
the  radiating  area  of  the  inside  of  the  trap,  the  net  heat  transfer  to  the 
window,  0  ,  may  be  written: 


where  d  is  the  diameter  of  the  window,  €  is  its  emissivity  for  long 
w  w 

wavelength  radiation,  Tw  is  the  window  temperature  and  T  is  the  mean 
trap  .temperature .  At  thermal  equilibrium,  Q^,  must  be  equal  to  the 
upward  heat  loss  from  the  window  (neglecting  the  small  amount  of  light 
absorption)  which  is  given  approximately  by: 


T 

w 


(269) 
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Using  these  two  relations 


T  - 
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*.  > 1/4 
1  * «.  j 


Hence,  the  upward  radiation  loss  is  given  by 

2 


ad 


%  " 
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which  is  much  smaller  than  the  heat  loss,  Q^,  for  a  simple  hole 
(Eq.  264)  since  is  small. 

The  total  upward  energy  loss  is  given  by  the  sum  of  the 

last  equation  and  the  reflection  loss,  rQ  ,  where  r  is  the  effective 

s 

reflectance  of  the  window  and  Qg  is  the  total  energy  input  to  the  hole. 
The  "breakeven"  point  occurs  when 


0.  +  ^  -  Qh 


or  where 


(272) 


(273) 


13.3.2  Laboratory  Test  Traps 

Two  fundamental  types  of  radiation  traps  have  been  assembled 
and  tested  during  the  current  Electro-Optical  Systems  program.  Examples  of  the 
wavelength  or  spectrally  selective  type  of  trap,  as  discussed  in  Section  13.3, 
are  the  cupric  oxide  on  aluminum  by  the  Unger  process,  the  deposition  of  black 
nickel  on  copper,  and  the  use  of  NESA  glass  with  a  blackened  absorber.  The 
multi- element  Fresnel  lens  fabricated  from  16  individual  vinyl  plastic  lenses 
(Bolsey  Research  and  Development  Corporation)  was  used  to  build  a  geometric 
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type  trap.  These  converging  lenses  focus  the  parallel  rays  of  incident  solar 
radiation  through  16  small  holes  in  an  aluminized  sheet  which  is  placed 
midway  between  the  multi-element  lens  and  the  black  absorber.  (The  radiation 
is  absorbed  by  the  collector  and  is  trapped  since  the  aluminized  sheet  has  a 
low  emissivity  and  reflects  collector  radiation  while  the  holes  are  so  small 
the  negligible  flux  passes  upward  through  them.) 

The  various  black  surfaces  produced  during  the  present 
program  were  laid  down  on  metal  plates,  four  inches  on  a  side.  A  test  jig 
was  designed  and  constructed  to  hold  these  black  surfaces  with  either  of  the 
radiation  trap  envelopes  previously  described.  When  using  either  the  NESA 
glass  cover,  or  the  multi-element  Fresnel  lens,  end  losses  were  minimized 
by  the  provision  of  highly  reflecting  side  walls  of  vacuum  deposited  aluminum 
on  glass.  The  radiation  trap  thus  consisted  of  either  a  selective  black  surface 
by  itself,  or  a  black  surface  discriminating  envelope,  and  the  reflecting  side 
walls.  The  trap  enclosure  is  thermally  insulated  from  the  metal  supporting 
frame  with  transite  spacers.  The  frame,  in  turn,  slides  onto  four  rods, 
attached  to  a  circular  metal  collar  that  fits  into  the  window  throat  of  the 
vacuum  test  chamber.  The  solar  flux  entering  the  test  chamber  window  falls 
normally  onto  the  surface  of  the  trap.  Equilibrium  temperature  measurements 
on  the  black  surfaces  were  affected  by  attaching  five  iron  constant 
thermocouples  to  the  surfaces.  The  back  surfaces  of  the  plates  have  a  low 
emissivity  and  provision  is  made  so  that  they  may  be  covered  with  a  layer  of 
insulation  if  desired.  Convective  losses  are  reduced  to  a  minimum  by  evacuating 
the  test  chamber  to  a  pressure  of  a  few  microns  of  mercury.  The  trap  structure 
is  shown  in  Fig.  51. 

13.3.2.1  Estimate  of  Thermal  Equilibrium  Conditions 

At  thermal  equilibrium,  the  temperature  of  a 
radiation  trap  may  be  calculated  by  solving  the  following  equation: 

4  4 

aA  q  *  P  +  K  +  A  €  or  T  -  A  OL  o  1  (274) 

sp^s  p  pbb 

where  a  is  the  absorptance  for  solar  radiation,  A  is  the  area  of  one  side 
s  P 

of  the  plate,  qg  is  the  solar  heat  input,  P  is  the  total  power  withdrawn, 
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K  represents  any  undesirable  thermal  conductivity  heat  loss,  c  is  the 
emissivity  (at  temperature,  T)  of  the  plate  and  a  is  the  Stef an* Boltzmann 
constant,  is  the  absorptance  of  the  plate  for  environmental  radiation 

whose  temperature  is  T^.  (Both  e  and  may  be  interpreted  as  the  sum  of 
the  contributions  from  the  front  and  back  surfaces  when  the  plate  is  not 
insulated  on  the  back  side.) 

The  heat  loss  due  to  the  thermal  conductivity 
of  the  air  may  be  calculated  knowing  the  geometry  of  the  experiment.  It  is, 
of  course,  a  function  of  T  as  well  as  the  environmental  temperature  T^.  Using 
this  Eq,  274  it  was  possible  to  correlate  the  experimental  work  done  in  the 
solar  test  facility  with  calculations  based  upon  the  measured  emissivity 
characteristics  of  several  radiation  traps.  (The  evaluation  of  equivalent 
absorptances  and  emissivities  requires  the  solution  of  integrals  involving 
Planck*s  function.) 

In  order  to  obtain  some  feeling  for  the 

temperatures  and  power  levels  involved,  Eq.  274  will  now  be  investigated 

under  the  conditions  likely  to  be  encountered  in  a  configuration.  The 

environmental  temperatures,  T^,  will  be  very  much  smaller  than  the  absorber 

temperature  T,  and  since  (X  =*  €,  the  last  term  in  the  above  equation  may  be 

b  4 

ignored  compared  with  €  a  T  .  Losses  due  to  thermal  conductivity,  e.g., 
through  support  structures,  can  be  made  sufficiently  small  so  that  K  may  be 
assumed  to  be  zero.  Thus,  it  is  assumed  that  the  absorbed  flux  is  partly 
transferred  as  thermal  energy  to  the  stream  and  partly  lost  by  reradiation, 
then 

4 

(X  A  q  -  P  +  A  €  or  T'  (275) 

s  p  s  p  x 


If  this  equation 
unit  area  of  the 


is  solved  for  the  power  transferred  to  the  stream  per 
absorbing  surface,  P/A^,  then 


P  4 

—  -  a  qs  -  e  a  T 

P  s 


(276) 
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If  a  collection  efficency,  is  defined  as  the  ratio  of  the  left-hand  side 

of  Eq.  276  to  q  ,  the  solar  input  per  unit  area  of  the  absorber,  then  the 

s 

ratio  of  the  collector  efficiency  to  the  absorptance  for  solar  radiation  can 
be  written: 


a 


a  q  A 
s  s  p 


1  - 


a 


a  T 


(277) 


These  last  two  equations  are  plotted  in 

2 

Fig.  52  for  qg  =  442  Btu/hr  ft  and  completely  specify  the  operation  of 

the  solar  absorbers.  Note  that  for  a  two-sided  collector,  €  represents  the 
sum  of  the  equivalent  emissivities  of  the  front  and  back  surfaces.  For  a 
case  where  the  collector  is  not  a  plane  surface  but  some  sort  of  airfoil, 

6  must  be  increased  to  account  for  the  fact  that  the  reradiating  area  is 
somewhat  larger  than  the  normal  area,  A^,  which  is  absorbing  the  flux. 

13.3.2.2  Experimental  Results 

The  unique  solar  energy  test  facility,  described 
in  Section  13.5  provides  a  means  by  which  solar  energy  collectors  can  be 
evaluated  throughout  a  range  of  flux  densities  and  air  pressures.  It  provides 
a  direct  means  of  supplying  the  following  information  for  various  collectors: 

1.  Stability  of  the  emissivity  and 
absorptance  characteristics  at 
high  temperatures. 

2.  Equilibrium  temperatures  at  various 
levels  of  operation. 

Stability  studies  have  been  carried  out  by 
adjusting  the  flux  input  (servo  controlled  screen)  to  yield  a  desired  temperature. 
This  value  of  flux  input  is  maintained  at  a  constant  level  despite  variations  in 
the  energy  incident  upon  the  searchlight  by  means  of  the  feedback  servomechanism 
controlling  the  screen.  The  temperature  of  the  sample  absorbers  is  monitored  as 
a  function  of.  time  and  displays  the  stability  of  the  surface  directly. 
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Equilibrium  is  achieved  when  the  difference 

betweien  the  absorbed  energy  and  the  energy  withdrawn  is  equal  to  the  energy 

lost;  the  equilibrium  temperature  for  various  energy  levels  of  operation 

can  be  determined  by  simply  varying  the  flux  density  on  the  absorber  surface. 

That  is  to  say,  it  is  not  necessary  to  attach  any  sort  of  heat  exchanger  to 

the  surface  to  withdraw  energy  in  order  to  study  equilibrium  temperatures  at 

various  levels  of  operation;  it  is  only  necessary  to  reduce  the  energy  input. 

The  solar  test  facility  has  been  in  operation 

fcr  some  time.  During  early  tests,  the  equilibrium  temperature  of  an  alum^ 

2 

inum  surface  with  a  cupric  oxide  coating  at  irradiance  of  .160  Btu/sec  ft 
was  found  to  be  350°F.  However,  it  was  found  that  the  flux  was  not  uniform 
over  the  surface  of  the  sample  because  of  optical  imperfections  in  the 
secondary  mirror  and  nonuniformity  of  the  flux  control  screen.  Furthermore, 
the  conductive  heat  losses  through  the  insulation  of  the  sample  holder  were 
larger  than  expected. 

Effort  is  now  being  exerted  to  improve  the 
system  by  redesigning  and  building  a  new  sample  holder.  In  the  interim,  it 
has  been  possible  to  take  meaningful  data  by: 

1.  Suspending  the  test  samples  with  nylon 
thread  in  order  to  reduce  conduction  losses. 

2,  Allowing  radiation  to  fall  directly  upon 
the  vacuum  test  chamber  without  concen¬ 
tration  in  the  mirror  system  and  without 
the  flux  control  screen  in  place  in  order 
that  the  flux  be  completely  uniform. 

The  last  factor  prevents  taking  measurements  at  other  than  one  level  of  operation; 

2 

at  about  0.080  Btu/sec  ft  net  input.  (This  would  have  the  same  net  effect  as  a 

2  2 
flux  input  of  0.119  Btu/sec  ft  and  a  useful  power  output  of  0.039  Btu/sec  ft  ). 

Both  the  black  nickel  on  copper  collector  and 

the  cupric  oxide  on  aluminum  collector  have  been  evaluated  in  the  solar  test 

facility.  The  latter  achieved  an  equilibrium  temperature  of  398°F  with  an 

erradiance  of  0.0845  Btu/ft^  sec,  while  the  black  nickel  surface  reached  452°F 

with  an  irradiance  of  0.0814  Btu/ft^  sec.  The  absorbing  surfaces  were  operated 

at  these  temperatures  for  thirty  minutes  with  no  apparent  detrimental  affects 

to  the  emissivity  characteristics  of  the  surface,  (Further  stability  tests  will 

be  made  when  the  problems  mentioned  earlier  have  been  corrected).  The  equilibrium 
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temperatures  measured  were  within  8  percent  of  the  temperatures  predicted  from 
Eq.  i 274  ,  the  errors  being  mainly  attributed  to  inaccurate  graphical  integration 
of  the  effective  emissivity  integrals  and  to  an  incomplete  knowledge  of  the 
diffuse  reflectance  of  the  surface. 

13  A  Interference  Filters 

The  optical  interference  filter  offers  perhaps  the  most  versatile 
available  method  for  controlling  the  spectral  radiation  properties  of  a  solid 
surface . 

Interference  filters  are  composed  of  one  or  more  layers  of  optical 
materials  which  have  varying  thicknesses  and  indices  of  refraction  and  which 
are  placed  upon  a  transparent  or  reflecting  surface  called  a  substrate.  The 
thicknesses  and  refractive  indices  of  the  films  are  chosen  in  such  a  way 
that  energy  of  some  wavelengths  will  be  largely  transmitted  or  reflected  and 
energy  of  other  wavelengths  will  be  cancelled  out.  Interference  filters 
involving  more  than  one  surface  layer  are  frequently  termed  multilayered  films. 

In  general,  the  filters  of  interest  for  solar  energy  collection  would  involve 
at  least  two  surface  films. 

By  arranging  the  layers  in  the  correct  manner,  a  filter  can  in 
principle  be  designed  with  a  specific  spectral  characteristic.  The  methods 
by  which  such  a  design  would  be  developed  are  the  subject  of  this  section. 

Although  interference  filters  with  very  good  spectral  characteristics 
have  been  developed,  other  types  of  absorbing  surfaces  offer  more  promise  for 
the  present  application.  Disadvantages  of  the  interference  filters  include 
their  probable  directional  sensitivity  and  the  difficulty  of  applying  them  on 
the  thin  and  fragile  surfaces  considered  for  this  application.  Therefore,  the 
emphasis  on  this  program  has  been  concentrated  on  other  types  of  selective 
surfaces.  However,  much  analytical  and  survey  work  has  been  accomplished  by 
the  Electro-Optical  Systems1  staff,  as  is  discussed  in  the  following  pages. 

13.4.1  Literature  Survey 

Technical  literature  has  been  reviewed  (Ref.  29  through 
Ref.  34)  which  disusses  both  the  theoretical  and  experimental  aspects  of 
Interference  filters.  The  literature  survey  indicates  that  interference 
filters,  particularly  of  the  reflection  type,  offer  promise  for  solar  energy 
collection  applications.  However,  it  is  apparent  also  that  much  remains  to 
be  established  concerning  the  general  design  theory  necessary  to  produce  an 
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arbitrary  desired  spectral  characteristic.  Methods  are  at  hand  for  analyzing 
a  particular  filter  with  given  film  thicknesses  and  optical  properties,  but 
the  direct  problem  of  designing  a  filter  with  prescribed  characteristics  is 
still  in  its  infancy. 

13.4,2  Theory  of  Thin  Film  Interference  Filters 

Filters  of  this  type  consist  of  a  series  of  metal  and 
dielectric  films  of  varying  refractive  indices  and  thicknesses  laid  down 
on  a  substrate  medium.  Their  principle  of  operation  is  thus: 

When  a  light  ray  traverses  these  layers,  the  new  rays 
formed  at  the  film  interfaces  have  controllable  phases  and  amplitudes  with 
respect  to  each  other.  By  arranging  the  thicknesses  and  indices  properly, 
the  new  waves  can  be  added  algebraically  in  a  wavelength  dependent  fashion 
to  achieve  the  desired  amount  of  cancellation  or  reinforcement. 

Interference  filters  are  classified  as  either  transmission 
or  reflection  types.  It  is  presently  felt  that  the  reflection  type  filter 
offers  the  greatest  promise  for  solar  energy  applications  for  the  following 
reasons : 

1.  No  transmission- type  filters  have  been  found  which 
possess  sufficiently  wide  band  pass  characteristics 
in  the  infrared. 

2.  Since  transmission  filters  are  directionally  sensitive, 
a  transmission  filter  would  be  partially  ineffective  in 
reflecting  the  diffuse  reradiation. 

The  reason  why  reflection  filters  are  promising  may  be  seen  by  considering 
the  general  construction  features  of  these  filters.  Briefly,  they  consist 
of  a  series  of  layers  arranged  in  the  following  manner: 

A  -  Substrate 

B  -  Highly  Reflecting  Metallic  Layer 

C,D-  Coatings  to  Produce  Low  Visible  Reflectance 
In  these  filters  the  combined  thickness  of  layers  C  and  D  is  of  the  order  of 
a  few  half  wavelengths  of  visible  light,  but  it  is  far  less  than  a  wavelength 
of  infrared  radiation.  Then  in  the  infrared  region  of  the  spectrum  the  phase 
retardation  of  the  reflected  rays  is  very  small  and  all  such  rays  add  in 
amplitude  producing  high  infrared  reflectance  and,  therefore,  low  infrared 
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emissivity.  At  present  there  are  two  ways  of  constructing  filters  of  this 

type.  These  are  depicted  in  the  table  below. 

Type  "A"  Type  "B" 

Layer  C  Semiconductor  ^o  Dielectric 

4  4 

Layer  D  Dielectric  ^o  Metal  (very  thin) 

4 

The  notation  means  that  the  layer  has  a  thickness  of  a  quarter  wavelength 

4  ^ 

for  some  visible  wavelength,  o. 

It  is  apparent  that  considerable  latitude  is  possible  in  the 
design  of  these  filters  since  and  the  materials  may  be  chosen  at  will. 

Also,  there  are  many  possible  variations  on  the  two  types  of  design.  For 
instance,  the  top  two  layers  may  be  repeated  several  times,  as  in  the  sequence 

ABCDCDCD  . ,  etc,  or  a  pair  of  dielectrics  may  be  put  on  top  giving  the 

sequence  ABCDE.  The  latter  should  give  an  additional  minimum  in  the  reflect¬ 
ivity  versus  wavelength  curve  and  lower  reflectivity  throughout  the  visible 
because  of  better  matching  of  the  top  film  to  free  space. 

The  design  of  a  filter  to  meet  the  required  specifications  is 
by  no  means  a  straight  forward  analytical  task;  but,  rather,  requires  a  long 
series  of  iterative  calculations  in  order  to  determine  how  the  best  use  may 
be  made  of  the  materials  which  can  be  deposited  as  films,. 

13.4.3  Three-Layer  Film 

There  are  several  basic  designs  which  deserve  consideration. 
However,  it  is  most  instructive  to  consider  the  simplest  of  these  designs.  A 
typical  design  and  reflection  characteristic  are  shown  in  Fig.  53. 

Among  the  parameters  which  may  be  varied  are: 

1.  Thickness  and  choice  of  metal  for  first  layer, 

2.  Refractive  index  of  dielectric, 

3.  Choice  of  metal  for  third  layer,  M9 

An  analytical  study  by  Hadley  and  Dennison  (Ref.  34)  yielded 
equations  for  the  performance  of  this  filter.  The  things  seen  to  be  generally 
true  were: 
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1. -  The  minima  in  reflectivity,  R,  could  be  brought 

down  to  zero  for  only  one  value  of  resistance 
(377  ohms  per  square)  of  the  top  layer, 

2.  This  resistance  is  measured  between  edges  of 
any  sized  square  of  deposited  film  and  is 
independent  of  other  parameters . 

3.  As  the  refractive  index  of  increases  from  unity, 
the  difference  between  X  and  X,  increases,  and 
the  local  maximum  in  R  increases. 

Thus  increasing  the  refractive  index  of  effectively  moves  the  cutoff 
wavelength  (Xc)  farther  into  the  infrared  when  Xq  is  kept  constant.  For 
filters  of  this  type  which  absorb  the  visible,  the  difference  (X^~Xq) 
only  about  0.6  4.  Thus,  since  the  difference  ought  to  be  more  like  2  4, 
a  higher  index  should  be  desirable.  However,  a  compromise  must  be  reached, 
since  a  higher  index  would  make  the  filter  less  absorbing  between  Xq  and  X^ . 
For  a  given  choice  of  metals,  then,  the  results  of  several  different  choice 
of  indices  must  be  compared,  and  the  best  of  these  results  for  each  metal 
must  be  chosen  as  the  optimum  filter. 
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14.  CONCLUSIONS  AND  RECOMMENDATIONS 

Without  considering  the  value  of  a  solar  powered  aircraft  of  con¬ 
figuration  and  performance  as  developed  in  this  study,  the  following 
conclusions  and  recommendations  are  offered  regarding  the  technical 
feasibility  of  such  a  device: 

a.  The  results  of  this  study  indicate  that  a  solar  powered 

aircraft  of  large  aspect  ratio,  1000-feet  wing  chord  length, 

flying  continuously  directly  beneath  the  sun  at  an  altitude 

2 

of  about  220,000  ft  with  wing  loading  of  0.42  lbs/ft  is 
feasible,  based  on  considerations  of  aerodynamics,  heat 
transfer,  and  energy  collection. 

b.  However,  before  undertaking  the  large-scale  experimental 
development  required  for  final  realization  of  the  solar 
powered  aircraft,  a  more  elaborate  design  and  performance 
analysis  is  necessary.  This  analysis  should  include  complete 
and  detailed  machine  performance  calculations,  as  well  as 
consideration  of  the  following  problems  which  were  not  in¬ 
cluded  in  the  scope  of  the  present  study:  structures,  control 
and  stability,  maneuverability,  launch  and  recovery,  dura¬ 
bility,  operation  non-normal  to  the  sun  direction,  cost, 
applications,  and  comparison  with  competitive  systems. 
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TABLE  I 


FLIGHT  CONDITIONS  FOR  TWENTY-FOUR  HOUR  EARTH  SATELLITE 
SOLAR  POWERED  AIRCRAFT 


H 

p 

T 

ft/1000 

i  Po 

°R 

190 

3.11 

-4 

10  4 

474 

200 

2.08 

457 

210 

1.36 

?4 
10  4 

432 

220 

8.71 

-5 

10 

410 

230 

5.44 

389 

240 

3.32 

367 

250 

1.96 

354 

260 

1.16 

354 

H 

= 

P/P 

0 

= 

T 

= 

M 

q 

= 

(v/V  )L 

= 

pr-tr 

Mil  >2 

m 

M 

q 

0/  V 

— 

lbs/f t^ 

ft"1 

1.43 

1.03 

3550 

1.47 

0.76 

2575 

1.50 

0.454 

1850 

1*54 

0.307 

1290 

1.59 

0.204 

895 

1.63 

0.131 

605 

1.66 

0.080 

384 

1.66 

0.047 

226 

PR 

tr 

(y/x>  )2 

0 

„  -1 

mm  Hg 

R 

ft 

0.738 

657 

1540 

0.522 

643 

1130 

0.356 

613 

806 

0.242 

594 

568 

0.163 

576 

395 

0.105 

552 

272 

0.065 

540 

173 

0.039 

540 

104 

Flight  altitude 

Ratio  of  ambient  pressure  to  sea  level  pressure 

Atmosphere  temperature 

Aircraft  fligh  Mach  number 

Flight  dynamic  pressure 

Aircraft  flight  Reynolds  number  per  foot 

Pressure  and  temperature,  respectively  for 
air  compressed  isentropically  from  free 
stream  conditions  to  M  *.3. 

Reynolds  number  per  foot  of  the  air  at  M  *.3 
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TABLE  II 


INPUT  DATA  FOR  PERFORMANCE  COMPUTATIONS 


*The  inadvertent  doubling  of  y^flm  in  this  case  doubles  the  friction  at  the  walls  for  the  internal  flow. 
**In  this  case,  the  friction  of  the  internal  flow  was  arbitrarily  halved. 


C-600-4 


1 


:nput  data  for  performance  computations 


Xc/y  X  /X-  T> 
f  p  r  ° 


Ere  E  E  E-  _  E 

ff  uu  11  uf  If  su 


Remarks 


.01142  0  I  .9150  .0200  .0300  .0300  .0200  .0200  .0308  y2/2m  doubled* 


.00285 


.01142 


.00285 


.01142 


used  1/2  f/2m** 


Not  completed 


>n  at  the  walls  for  the  internal  fiow. 
Lived. 


ft 

4 


TABLE  III 


SUMMARY 

OF  RESULTS  OF  PERFORMANCE 

COMPUTATIONS 

X 

Excess* 

.  ** 
w/x 

c 

Thrus  t 

c 

Case 

Feet 

Coefficient 

lbs/f 

I 

200.7 

-0.004939 

— 

II 

800 

-0.001651 

II 

1000 

-0.001600 

-- 

III 

300 

-0.008370 

III 

1000 

-0.002060 

III 

1400 

-0.001455 

III 

1800 

-0.001676 

-- 

IV 

400 

-0.003127 

-- 

V 

1200 

+0.001078 

.045408 

V 

1400 

+0.000891 

.041279 

V 

1430 

+0.000346 

♦025733 

VI 

1000 

+0.000230 

.020972 

VI 

1400 

+0.000139 

.016283 

VII 

1800 

-0.000299 

— 

VIII 

800 

+0.000214 

.022236 

IX 

Not  Completed 

X 

800 

+0.000385 

.027125 

XI 

800 

-0.001804 

XI 

1400 

+0.000254 

.022023 

XI 

1800 

+0.000850 

.041797 

XII 

1400 

+0.000489 

.030580 

XIII 

1800 

-0.000166 

-- 

XIV 

500 

+0.002689 

.032556 

*  Excess  Thrust  Coefficient  -  -  c  -  C 

qx  Do  ^Df 

o  c 

**  Wing  Loading  ■  w/x^  lbs/fc 
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FIG.  1  SOLAR  POWERED  AIRCRAFT  CONFIGURATION  -  WING  CROSS 
SECTION 


Note:  The  trends  shown  by  these  curves  are  correct.  The  wing  loadings 
shown  were  calculated  by  the  simplified  method  of  TN  59-184  and 


1 

FIG.  2  WING  LOADING  VS  ALTITUDE  FOR  VARIOUS  VALUES  OF  CHORD  j 

LENGTH  (FROM  TN  59-184) 
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i 


i 


23 


NOTATION  FOR  FIGURES  4  THROUGH  15 


F 

f2 

M0 


P* 

s 


s2 


thrust  per  unit  span,  lbs/ft 

~  x  (energy  flux  transferred  to  the  internal  stream) 
qs 

internal  stream  Mach  Number 

P 

—  =  internal  stream  stagnation  pressure  ratio 

s2 

internal  stream  stagnation  pressure 

internal  stream  stagnation  pressure  at  Station  2. 


rr*  — 

_f  = 

f|5 

heat  transfer  foil 

temperature 

ratio 

I 

1 

s2 

II 

H 

Ti  = 

lower  wing 

surface 

temperature 

ratio 

Ts2 

T* 

T 

u  - 

upper  wing 

surface 

temperature 

ratio 

u 

x*  =» 

s 

T  o 
ijSc 

s  = 

internal  stream  stagnation  temperature  ratio 

T 

s2 

Tg  “  internal  stream  stagnation  temperature 

X  -  internal  stream  stagnation  temperature  at  Station  2. 

S2 


T^,T^,Tu  =  temperatures  of  heat  transfer  foils,  upper  and  lower  wing 
surfaces,  respectively. 
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DETAILS  OF  PERFORMANCE 
COMPUTATION  RESULTS,  CASE  I 


FIG.  5 

DETAILS  OF  PERFORMANCE 
COMPUTATION  RESULTS,  CASE  II 


too 


ooo 


1200 


♦400 


MOC 


OF  PERFORMANCE  COMPUTATION  RESULTS, 
I 
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NOTATION  FOR  FIGURES  16  THROUGH  19 


insolation, 


2 

ft .lbs/sec .ft . 


Q 


heat  transfer  rate  per  unit  span  to  the  internal  stream, 
ft .lbs ./sec .ft . 

“x  (solar  energy  flux  absorbed  by  the  collector- transfer 
%  section.  ) 


C  =  heat  flux  by  convection  from  the  wing  surfaces  to  the  external 
ex~  stream,  ft .lbs/sec-ft^ 

C.  -  heat  flux  by  convection  from  the  collector- transfer  foils  to 
the  internal  stream,,  ft  .lbs/sec-  ftz 

C.  «  heat  flux  by  convection  from  the  collector- transfer  wall 
3.  n  t  *  o 

surfaces  to  the  internal  stream,  ft .Lb/sec-f t . 


x  =  chordal  distance  from  wing  leading  edge,  ft. 

R  m  radiation  flux  lost  from  upper  and  lower  wing  surfaces,  f t . lbs/sec- ft 
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100  200  300  400  500  600  700 


FIG.  16  DETAILS  OF  PERFORMANCE  COMPUTATION  RESULTS, 
CASE  XIV 


^  dx 


FIG.  17  ENERGY  FLUX  DETAILS,  CASE  VIII 
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S  *  7*  %  s  s  <. 


FIG.  18  ENERGY  FLUX  DETAILS,  CASE  XI 


_L^£.  =  -n  _ 

a*  *«  q,  a,  q,  q. 


0  200  400  600  900  1000  1200  1400  1600  »00 


X 


FIG.  19  ENERGY  FLUX  DETAILS,  CASE  XII 


i 
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FIG.  20  INLET  DIFFUSER  CONFIGURATION 

q*  \ 


FIG  21  ARRANGEMENT  OF  RADIATION  COLLECTOR-HEAT  TRANSFER  FOILS 
IN  THE  COLLECTOR-TRANSFER  SECTION 
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qu 


-UPPER  WING  SURFACE 


FIG.  22  ELEMENTARY  RADIATION  COLLECTOR  CONSISTING  OF  INFINITE 
PARALLEL  SURFACES 


FIG.  23  RADIANT  ENERGY  EXCHANGE  BETWEEN  INFINITE  PARALLEL  SURFACES 

I" 


FIG.  24  "GREEN  HOUSE"  EFFECT  BETWEEN  INFINITE  PARALLEL  SURFACES 
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TEMPERATURE  -  dtflrot ,  R  ANKINE  TEMPERATURE  -  d««rs«s ,  RANK1NE 


FIG.  27 

AIR,  PRANDTL  NUMBER 


a)  Relation  of  absorber  characteristics  to  rodiation  spectra 


FIG.  28 

ABSORPTION-EMISSION 
CHARACTERISTICS  OF  SOLAR 
ENERGY  ABSORBING  SURFACES 
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FIG.  29 


\  \ 

\  -  TRAPPING  EFFICIENCY  (f)  VS 

\  CUTOFF  WAVELENGTH  (\  )  FOR 

\  _  AN  IDEAL  COLLECTOR  C 

J _ I _ I _ L _ 

1.5  2.0  2.5  3.0  3.5 

Xc-microns 


FIG.  30  SPRAYED  COATINGS  ARE  EVENLY  APPLIED 
TO  TEST  SAMPLES  BY  USE  OF  CAM- 
OSCILLATED  NOZZLE  AND  TRAVELING 
HOT-PLATE 


REFLECTANCE -p«rcenl 


FIG.  33 

TRANSMITTANCE  OF  NESA  GLASS, 
0.25^-1411  (SAMPLE  NO.  P-6) 
10.5  ohms/sq. 


FIG.  34 

TRANSMISSION  CHARACTERISTIC 
OF  GOLD  SMOKE  FILTERS 
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TRANSMITTANCE-  percent  ABSORPTIVITY- p*re*nt 


FIG.  35  OPTIMUM  EMISSIVITY  VS  ABSORPTIVITY 
REFERENCE 


FIG.  36 

TRANSMITTANCE  OF  PYREX  AND 
NESA  GLASS 


TRANSMITTANCE  OF  PYREX  AND  NESA  GLASS,0.25,»-2.!* 
fitunpl*  No.  P-6,tMckMW^ 
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REFLECTANCE-  p«rc«nt 


FIG.  37 


CUPRIC  OXIDE  ON  ALUMINUM 
BY  THE  UNGER  PROCESS 


FIG.  38 

CUPRIC  OXIDE  ON  ALUMINUM 
PREPARED  BY  THOMAS  A.  UNGER 
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FIG.  39  REFLECTANCE  OF  UNCOATED  ALUMINUM 


FIG.  40 

SPECTRAL  RESPONSE  OF  CUPRIC  OXIDE 
ON  COPPER  BY  THE  "EBONOL"  PROCESS 
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0  I  2  3  A  5  6  7  8  9 

WAVE  LENGTH -microns 


FIG.  41  BLACK  NICKEL- IRON  OXIDE  ON  STEEL  BY  THE 
"NICKEL  PENETRATE"  PROCESS 


FIG.  42 

BLACK  NICKEL  ON  COPPER 
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REFLECTANCE-  prrctnt 


f 

j 


FIG*  43 

BLACK  NICKEL  ON  COPPER  SHOWING  THE 
EXTREMES  IN  SPECTRAL  CHARACTERISTIC 
0. 35^-13^1 


FIG.  44 

BLACK  NICKEL  ON  NICKEL  AND 
NICKEL-PLATED  COPPER 


4 


4< 


FIG.  45 

PLATINUM  BLACK  ON  PLATINUM 


FIG.  46 

CARBON  BLACK  ON  ALUMINUM 
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FIG.  47 


REFLECTANCE  OF  NESA  AND 
CORNING  GLASS 


CORNING - ll.Ooi.ms/sq. 

NESA  ^9 - 6.6  ohms/sq 

NESA  *28 .  20.0  ohms/sq. 

NESA  *27 -  14.7  ohms/sq. 


5  6  7  8  9  10  1 1  12  13  14 


WAVELENGTH -microns 


FIG.  48 

REFLECTANCE  OF  NESA  GLASS 


4 


REFLECTANCE -P«rc«nt  (  TRANSMITTANCE -percent 


FIG.  49 

TRANSMITTANCE  OF  NESA  GLASS 
THICKNESS  1/8" 


FIG.  50  TRANSPARENT  SURFACES 


FIG.  51  FRESNEL  LENS  AND  NESA-GLASS  RADIATION  TRAPS  FIT  INTO 
TEST  FRAME  FOR  MOUNTING  IN  SOLAR  TEST  FACILITY 


FIG.  52 

SOLAR  COLLECTOR  CHARACTERISTIC 
FOR  SPACE  APPLICATION 


(a)  Typical  Design 


(b) Typical  Reflection  Characteristic 


FIG.  53  THREE -LAYER  INTERFERENCE  FILTER 
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